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abstract 


In this thesis entitled A UNIFIED APPROACH FOR 
MODELLING BJT BASED ON GUMMBL POON MODEL a set of 
parameters related by a sert of model equations has been 
postulated which specify transistor electrical behaviours 
for all regions of its operation 

The device model chosen is based on Gummel-^oon 
charge control model Extensions/modifications are made 
wherever necessary Help is taken from some recent 
techniques reported for modelling nonlinear effects that 
occurs insich a transistor In the remaining cases we 
have made our own derivations and also that r quired to 
accommodate various nonlinear effects occurring simul- 
taneously 

The parameters are extracted using nonlinear 
least square curve fitting algorithm suggested by 
Levehberg-^arquardt The data points are obtained from 
the experimentally observed char aot eristic curves of the 
transistor To take measurement a precision curve 
tracer is built as a part of the project work 

The parameters chosen are physical parameters 
and their relationship with device material and structural 
Parameters are indicated These relationships can be 
utilised for checking the validity of the model 

Finally con^arision between the proposed model and 
existing SPICE model is made 
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INTRODUCTION 

1 1 THE CONCEPT OF MODBDliING 

ihe tdim model has a number o£ overlapping meanings 
which can easily create oonf\i3lon The most common meaning 
encountered is the refeirence to the simulation of the pl^si- 
cal appearance of an object generally on a different physi- 
cal scale A more relevant meaning for our puroose la the 
duplication of its physical appearance one then can cons- 
truct a set of equation e a mathematical model) to 
portray the internal behaviour as well as that with the surr- 
oundings A still more nDbust meaning has to do with the 
as if analysis of tlie system By noting some external 
behaviour using perturbation one replaces the system (a 
black box) by a set of interacting elements with known 
behaviour This conceptua3i.«ed (as opposed to constructed) 
model presents circumstances with knoVn principles which 
acts as a generator of p^^inciples in the theorist thinking 
The laws or operating principles in the model postulated 
are assumed to hold good for tha system The laws for the 
system being so modelled constitutes the theory It is also 
possible to go further and assume that the physical consti- 
tution of the unknown ci'^curostances requiring a theory is 
similar to that of the model (vis the hybrid model of a 
transistor) ihe use of either hypothetical syeteims or 
analogous physical systems allows for eidh conjecture which 
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in i:um allows more powerful theories tp emerge The analogy 
of the physical model builds the theory for the model or if 
the theorist recognizes that he has only a partial analogy 
in his model the model at least contributes some axioms to 
his theory (e g the example cited above) 

\ 2 AIM OF the work 

It is also poss'^ble to build theories as opposed to 
conceptual or physical models and still obtain some end 
result But the advantage of models is that they usually 
make theoretical development easier Though the simulation 
based on Gunntel Poon model in its early days has incorporated 
both physical and conceptual theories with more and more 
understanding of the behaviour of a transistor available 
nowadays it is possible to reformulate the model fully based 
on physical theories The aim of this work is towards this 
line Gf course prior validation of the theory and the 
verification of the software package are mandatory before one 
can ude it to study or simulate larger systems There are 
varioiis levels of modelling we have more to say about this 
in Chapter 4 

1 3 ORGANISATION OP THE THESIS 

The structure of this thesis is as follows In 
Chapter 2 we have started with the Gummel Poon model the 
focus of the work Later on various auxiliary effects that 
significantly affect the functioning of the transistor are 
discussed some more parameters are also included to conaide| 
some new effects which are not taken care of in the Gtunmel Po<jk 



ntodel initially proposed we have also replaced some •£ 
the old theories by new theories recently developed by others 
or by us In Chapter 3 we have gone to develop a curve tracer 
to take measurements srequired for the evaluation of the model 
parameters In Chapter 4 in the introductory part we have 
reiterated the concept of modelling with its relevance to 
our work In Section 4 2 we have made a systematic discus- 
sion about how we can find out the parameters from some 
typical measurements whereas in the very next section we 
have hinted about the salient points in out proposed model 
Finally we have indicated the con^tibiUty of ours with 
the SPICG model and merits demerits between them 
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CHAPTER 2 

GUMMBIj POON MODBIi 


2 1 INTRODUCTION 

Any model on bipolar junction transistors must reflect 
the fact that a transistor consists of two p-n junctiens The 
simplest yet novel version of it is the Ebers-Moll model [6 ] 
where the two junctions are represented by ideal exponential 
functions relating current and voltages This modal erd^odies 
superposition that is that the transported cun^snt from 
emitter to the collector can be expressed as a sum of two 
independent conqponents resulting frrm the biasing of the 
junctions or the junctions act independently For real 
transistors the violation of superposition principle is eaallj 
observed for example the Early' effect [37] denotes the 
dependence of low frequency output conductance on bias For 
a transistor operating in its active region (Figure 2 1(d)) 
the collector current is seen to increase slowly with 
collector-emitter voltage v^^ This is due to base width 
modulation by the GOllector'4)aae bias voltage Ebers-Moll a 
model is unable ts include this and other second order effect^ 
like high injection effects bias dependent current gain 
An improvement in this respect is achieved by the Gtomme 1-Poor 
Model [7 ] based on charge control theory an extension on 
Mbll-Ross version [s] Uptil now Gummel-Pooa modal is widel* 
accepted for designing transistors as it is capable of 
describing the second order effects (Just mentioned above) | 
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more or less accurately However we will show In our work 

this model is still capable o£ handling some other important 

effects like conductivity modulation in the base and in the 

low-doped collector regions emitter crowding etc 

Since our work is basically based on this model 

we will first describe it in detail side by side we will 

point out the approximations taken at different stages and 

their validity for variotis region of operation of transistors 

thereafter we will passon to derivation of base-transit time 

as a function of current minority carrier chacges in the 

base region and finally emitter crowding effect 

2 2 1 BASIC DERIVATION OP G-P MODEb EQUATION 

At any point in a semiconductor electron and hole 

current densities and i are related to the carrier con- 

•'n •'p 

centrations by the expression 

(2 1 
•t 

and jp = q IJ’P P B - q Dp V p 

Multiplying the first one by and the second one by ji^^n. 
and then substracting so as to eliminate B we end up with 

l*p P Jn “ ‘‘n " ^ 

where Einstien s relation D ® p.kT/q for either type of 
oarriers is assumed Ihis expression is quite general and 
in one dimensional form it reduces to 

P j„ - n Jp = li^, Ujp W §5 (np-) (2 2i 
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Let \is apply this result in the base region of a 
transistor Gummel has assuraed the base currant to be negli 
gible i e unit emitter efficiency and no re^omoination 
with this assumption * 0 and » const » J„ (called the 

''n CO 

transported currant) and equation (2 2a) becomes 

P^CC “ ^Vl ix 3 


Integrating this expression from Xg the edge of the base 
depletion layer on the emitter end to x ® x_ the edge of 

cs 

the depletion region on the collector end we get 


pi^ 

/ C( 

Xg 


kTCrjgPg - n^p^) 


(2 4 


where suffix B and *€ is related to quantities at the 
emitter and collector Junctions respectively If we now use 
the junction law relating charge densities with junction 
voltages as 


2 ^bb'^T 


and 




„2 /bc/^T 


(2 


we will arrive at 


j 


CO 


kT n 


2 e 




e 


^BC^T 




/ p dx/yL 


(2 k 


which is the Guinmel equation Ihe prime sign is used with 
the junction voltages so as to differentiate from measured 
voltages If we neglect position dependant mobility the 
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denominator of (2 6) can be replaced by the where Qg is 
the majority carrier sheet charge density in the base region 
On further defining the following relationship 


'cc 


- 


% 


2 2 


^3 


Qj 


'e 


(2 7) 


where is the emitter area; from (2 6) the expression for 
current density can then be changed over to an expression for 
linking current 


’cc 


- T 


V /v 

. a 


(2 8 ) 


The same expressicm is derived by Gummel [31] in a somev^iat 
different fashion Here stands for the doped charge in 
the quasi neutral active base region when no h is applied 
and the quantity Ig is called the saturation current of 
the transistor 

Prom the modelling point of view two important 
featuxres should be noted which show the novelty of the G-P 
model Firstly the rapidly varying exponential factors 

VnipArp ^ ^ ^ 

e ^ and e appear explicitly in the n\imarator of 

equation (2 8) These factors do not depend on any model 
parameter (whereas in the B-M model the exponential terms 
depend on emission coefficients which are experimentally 
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detexnti.ned model parameters) Thus vjhen the junction voltages 
are known they can be evaluated accurately The remaining 
variable in (2 S') Is the total mobile charg*?^ n the base 
region which is more slowly varying function of bias than 
the exponentials in the numerator This leads to second 
important featuj^e of this new charge control relation model- 
ling of the bias dependence of the There is a great deal 

of flexibility involved in the modelling of 9b Trade offs 
are pcresible between complexity and accuracy In their ori- 
ginal paper [7 ] Gummel-Poon writes as the sum of the doped 
charge in the effective base region and the excess majority 
carrier charge The latter^quated with injected minority 
carrier chasrge by assuming quasi neutrality condition 
(n^(x) * Pj(3c)) It is further related to forward and reverse 
ourrents via charge control theory The break up is as 
follows 

* q Ag / (Nj 5 ^{x) + Ng(x))dx 

» q A / N^(x)dx + / ng(x)dx 

“ ^ «be ® Vf Vr 

The first term in the expansion represents the base 
charge at aero bias 0^^ represent the increase/decrease 

of base charge because of s^aduction/extension of depletion 
region in the base under the application of bias voltage in 
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BE and BC junctions respectively Their values are positive 
for positive bias and negative for negative bias in npn 
transistor The fou3Cth and fifth terms rep* ^ uhe Injected 
charge from the emitter and collector junction respectively 
In deducing this the so-called charge control principle is 
used which says that a stored chasrge Q is proportional to the 
related current I (Q « ti) where t has the physical signifi- 
cance of the transit time Thus and represent the 

forward and reverse transit times where I- and denote the 

X r 

forward and reverse currents In the fourth term the factor 

B is Introduced in order bo tackle base push out effect we 

will elaborate this point later on 

2 2 2 Approach in Modelling 

Let us now look in detail the intrinsic operation of 

a transistor In Figure 2 1 is shown a schematic side view 

of a discrete n-p-n transistor along with the flow of electronj 

and nole currents The solid line represents the flow of 

electrons and the dashed line represents the flow of holes 

The superscript e and n are to distinguish particle current 

symbols The collector current I consists of four component 

c 

“ ^cc - ^ol - ^c2 

where I is the dominated collector current through the tran 
sistor and its bias dependence is given in (2 8) It repre- 
sents the injection of electrons from the emitter into the 
base ^Ich traverse the base region (by diffusion and partly 
by drift when base grading is present) and ultimately collect 



LnlN(x)| 



Fig 21 Schematic cross section of a transistor (not to 
the scale) showing various current components 



Fig 2 2 Typical impurity profile for an epitaxial transistor 
showing charges associated with various regions 
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at the collector and small components due 

to recontoination in the collector-'base junction depletion 
region and in the bulk is current comm*=*nt generated 

by impact ionization Electrons generated by impact ionisation 
flow into the collector and become an additional current 
component to the collector current For a transistor biased 
in the active region is generally much larger than the 
other three The current is in fact important only at 
high CB junction reverse voltages Since we are not inter- 
ested in high voltage operation we will neglect its presence 
all throughout Modelling impact ionisation is a separate 

issue and we refer to [8 9] 

I comprises of two parts and already intro- 
duced in (2 9) represents a forward injected current due 

to electrons injected from emitter to the basa (for an npn 
transistor) ij, is the inverse electron current due to 
electron injection from collector to the base They can be 
expressed as 


I . (e ^ - 1) 


- 1) 


T s« 

r % 


(2 ll! 


Implicit is the assumptiorr of reciprocity theorem which is 
valid urs^er reasonable conditions Thus 


I 


GC 



I 


r 


(2 
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The total recombination rate is modelled phenomenologically 
and consists of five components two ideal and three non- 
ideal In absence of this also aquals th bc.se current 
Ijj given by 


\ “ ^el ^be2 ^e3 ’** ^bcl ^c2 

represents the holes coming in from the base contact 
that are injected into the emitter region These holes 
recombine in the bulk emitter region as shown in Figure 2 1 
Thus has the ideal exponential dependence on of 

the form 


I 

^el ^ " 1) 1^) 

’ F 


pp is called forward amplification factor This defi- 
nition is justified at low injection active region whence 
= %>o if one neglects the recombination generation 
current then this factor stands for the ra i^o I 

is related to device parameter by 

_ *e ** **8 . 

Pp 'itri^eoWv^/I^T ■ 


here Ig Ng stands for diffusion length depth and doping 

density in the emitter region respectively 

^be2 represents holes injected from the base contact 
which recombine In the emitter junction it can be modelled 


as 




lqji(S 




- 1 ) 


(2 15 
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Here is called space charge limited saturation current 
and Hg is a constant whose value lies between l and 2 An 
approximate expression for igg is 

^SE " « Ae ”i ®max <2 15a 


where x is the life time of electron or hole in the depletion 
region and is given by 


max 


2q - Vb^) 1/2 

*• e (N^iXg) + Ng) J 


(2 ISb 


vdiere » built in voltage N^CXg) is the doping level at 

the edge of the depletion region in the base Needless to say 
these quantities are related to BB junction 

The third component of (2 13) takes care of recombi<* 
lation in the bulk region of the base Using charge control 
theory it can be modelled as 



^3 - 


(2 16 


where is the recombination life time in the base region 
Using (2 9) we see that ^ In 

modem days planar transistor having narrow base width this 
ratio is quite small and one may neglect it as well 

The last two contributors to consideired in (2 13) 
comes from holes Injected from the base contact which recom- 
bine in the collector epilayer and Junction regions As in 
(2 14) and (2 IS) we can similarly model them 8 

<K 
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and 


bc2 


SC 




- n 


C2 17a) 


Here is the reverse amplification factor Ig^ is reverse 
saturation leakage current They can be attributed similar 
theoretical expressions as in (2 14a) and (2 15a) with the 
only difference is that the suffixes related to EB junction 
have to be replaced by that of collector base junction 

It is to be noted that in equatl n (2 13) we have 
not considered contribution from overlap diode comprising of 
inactive area of the CB junction This is important when 
CB junction is forward biased and we will consider it later or 
(Sec 2 7) .Equations (2 8) (2 9) and (2 13) sumnarize the basic 

equations that describe the intrinsic part of the transistor 
A similar description is possible for pnp transistor also 
2 3 NONLINEAR EFFECTS 

The nonlinearity of a translator as it appears in the 
basic E-M model via the exponential fom of the junction vol- 
tages in a way does not violate superposition as the cursrents 
derived from them are independent of the other But in a real 
transistor many secondary effects arise Ihey are listed as 
follows 1) Early effect 2) conductivity modulation in the 
base and collector epitaxial region 3) base push out effect 
4) emitter or edge crowding 5) degeneracy in the emitter 
region 6) voltage dependence of junction capacitances 
7) impact ionization 8) surface leakage currents and para- 
si tics Apart from these some physical phenomena like band 
gap narrowing# mobility behaviour of minority carriers gradlr^^ 
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of doping concentration in the emitter and base region varia- 
tion of transit time with current density etc to some extent 
impress on the functioning of a transistor some of which 
we will also take care of In the above list we will set 
aside the seventh one as it is important only at very high 
reverse junction voltage and we rather refer to [8 9 ] Let us 

discuss these effects critically in the light of Q-P model and 
modifications are suggested wherever necessary 
2 3 1 Earlys Effect 

The Early effect [37 17 20 ] arises because of the 

change in the effective base width as the collector-base Junc- 
tion width changes This effective decrease of base width (for 
reverse bias) rediaces the recombination rate in the base region 
(and hence in^rovement in transport factor) and more import- 
antly enhances the diffusion current because of sharper 
minority carrier gradient in the base region Prom a charge 
control point of view changes in Q_ the charge capacitlvely 
stored in the collector base junction cause changes in the 
collector current This give rise to an output conductance 
approximately given by 


Q 


out 


I?c 




<2 18) 


where is the collector current if there were no Early effec 
C CB junction capacitance 

The conventional Junction capacitance is related to 
junction voltage through an expression containing three para- 
meters 
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C 


oe 


(1 - vAbi) 


m 


(2 19) 


The parameters Sire junction b\illt-in voltage (typically 

^ 7 V for Si) the grading coefficient m the constant in the 
numerator which can be related to the zero bias capacitance 
But this expression suffers from numerical instability as v 
approaches C goes to Infinity In a real translator 

of course a finite amount of charge is stored for all bias 
conditions and the derivatives of chargee with respect to 
junction voltage is finite in order to avert this problem 
Guinmel has introduced a foujrth parameter in the denominator 
of (2 19) They have defined a function [39^ 


f(v P) 


P- 


(1 + P4> 


P. 



<2 IQa) 


where the vector p denotes the four parameters p 2 

p^ in brevity The normalised emitter and collector charges 

are then given by 



f {. ■ p ) 

(2 

20) 

%C “ 

^BC 

(2 

21) 


Though our software is based on this model concept a few 
points should be brought in before we pass on to the next 
topics Because of the complex nature of the function (3 19) 
the parameters are not readily amenable to numerical evaluation 
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Further in deriving these expressions the nature of the 
doping profile at the Junction is not considered and the 
expression for CB and eb junction comes undur th same token 
Since they are purely model parameters they are not as such 
directly related to device design parameters Hence evalua- 
tion of them does not provide any significant feedback to the 
designer So to be realistic we prefer to present here 
two different approaches The first one is due to us and the 
last one is somewhat close to the prescription [20] 

In Figure 2 2 is shown the doping profile and charge 
associated with various regions commensurate with our defini- 
tion The points and Xj^ specify the metallurgical 
junctions The points and represents the edges of 
depletion layers in the base when no bias is applied The 
points Xg and shows the same under a typical operating 
situation of the transistor So accojcding to our definition 

^co 

= Ag q / N^(x)dX 
*Eo = 0 

X 

“ Ag q / IJj^(x)<ix (2 22) 


Let U8 assume doping profile in the base region as 




A “CyX 

w N- e 
A 


for 0 i X i XjQ 


(2 22 ^ 

I 


-c lx I 
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where Cp and Cj^ can be called forward and reverse grading 
coefficient in the base They are both positive and their 
value can be approximated by 


S' 

s 


In 


A 


In N. 
A 

*jE 


(2 22b) 


we now can get closed form expression for and 0^^ 

Prom (2 22) and (2 22a) 


So 

Sc 

Se 


q ^ 




q Ag(s ^ ^ - e 


q A_j(l - e 


'r*S* 


(2 22c 


This is the situation in forward injection For reverse 
injection with bb in reverse biased only the expression for 
ought to be changed by simply replacing by In 

normalised form 


Sc 


Sc 

So 


- 1) 





3^ 

1 


■ e 


(2 22cr 


To relate these quantities with junction voltages is a fairly 
involved task (sea Appendix II) But this can be easily 
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computed from the constant collector epl -profile If we 
suppose an abrtipt space charge region with Identical amount 
of opposite charges on the both sides of the c Hector base 
junction region then 




(2 23) 


vrtiere is constant epi layer concentration w and v/^ are 
defined In Figure 2 2 

Assuming the base is much more heavily doped than 
the collector so that almost the whole voltage drop occurs 
across the collector depletion region then w can be approxi- 
mated as 


w 


L gj J 


and 


w. 


!■ ' J 


ep 


ep 


(2 24) 


Here is the collector junction built-in voltage one 

can introduce a model parameter m (say) instead of using half 
power factor in (2 24) to get better approximation Combi- 
ning (2 23) and (2 24) we get the ratio 


^30 



(2 25) 


as a function of with 

V « Q? /2<I ( 2 2< 

When approaches the Vp^p value the i^ole of the base 
region is depleted so we call this voltage the forward base 
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punch through voltage 

The above formulation Is valid so long as depletion 
region does not reach the heavily doped <n ) substrate When 
this happens we define the corresponding CB junction voltage 
as the epl layer reach through voltage Designating it by 
we can show 

rep 


V 


rep 


JL 

2e 




V 

'^BIC 


(2 27) 


For voltages greater than depletion width in the 

collector region practically stay® at w^p and junction capa- 
citance is then almost constant The stored charge in the 
space charge layer then linearly Increaseswith and can 
be fcmmlatad by 






^ ^ep V 


(2 26) 
(2 29) 


The ratio In this case would be (from (2 27) and (2 28)) 




* * ^rep^ 




V 

rep PlP 


(2 30 ) 


\^lch simplifies to 


Sbc 






PEF 


PBF 


> " ^PBP'^PIP ^ ^ 

(2 31) 



\diere 
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If approaches the base region becomes fully depleted 

Co PBF 

So can be called as the extrinsic forward base punch 

PEP 

through voltage 

Under normal operation the base emitter junction is 

forward biassed such that t* ^bic *" ^BIC ^ 

Equations (2 25) and (2 31) then reforms to 




vTv ' Aj •• Vy A/ 


(2 33) 


for v^g 


< V 


rep 





*** ^BIC 


[1 + 
^PEP 


(1 - 


^CE > 
^PEP 


(1 - 


^PEF ^ 

v"'" * 

PIP 



(2 34) 


Thus only two parameters and Vp^p (or are suffi- 

cient to provide fairly accurate modelling of Early effect 
\^er 0 as G-P requires four parameters to be evaluated 


The chaacge can be considered in a simple way 
Assuming one sided abrupt junction the depletion width will 
be approximated by 


X 


1 


^ YbIE 

«t<NAE> 


) 1/2 
-] 


(2 35) 


Hhare <NAB> is the average concentration of the inpurities in 
the base near the emitter defined by 


/ N^(x)dx 

<NAE> « 36J 

^ ^BO 

Of course <NAE> shows strong depends on However as 

is clear from the Figure 2 2 this average concentration does 
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not vary significant in typical reverse operating inode (one 
can replace it by na( 0 ) ) but shows strong dependence in the 
forward acting mode because of the reverse g'^ading effect 
close to But in the forward direction the variation of 

Vgg itself is not significant so this approximation cannot 
in any way impair the result so we can use an expression 
similar to (2 25) also for the ratio 




(v + V y/V - V yF Ar 


where 


PIR 


IqT-<NAS^ 


(2 37> 

(2 37a) 


we note from (2 37) that for « 1 Vp^j^ is the base 

emitter voltage at which the base is fully depleted So as 
before let this parameter be defined as rav-erse base punch 
through voltage 

2 3 2 Mobility Variation in the Base Region 

While going from ('’ 6) to <2 7) the mobility within 
the integral in (2 6) was replaced by an average value out- 
side the integral xn his original derivation [31} Gunsmel has 
shown that the variation of mobXlity gives rise to an effective 
widening of the base (or an equivalent increase in Gummel 
nunber) by an amount fii 0 02 Pm For modem days 

transistor this is significant and cannot be simply ignored 

The reason for this can be seen physically as follows 
Because the net in^urity concentration is not constant across 
the base region the peak of the profile occurs close to BE 
Junction where holes have got greater md3ility than those at 
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the falling edge of the pjroflle close to the CB Junction 
Thus as majority carriers are depleted from the base as a 
result of an Increase in it is the carriers with lower 

mobility weighting that are lost This results in a smaller 
fractional decrease of effective base Gunvmel number and hence 
less pronounced dependence of over than vrould be pre- 
dicted from charge control approach 

The single analytical es^ression of mc^ility in a 
semiconductor in its sin^lest form can be written as 











(2 38) 


2 

where - 1400 cm /V sec is the mobility of intrinsic 

silicon a = a 72 “ 8 5 El6/cm^ and « 7 5 x 10^ V/cm 

(critical field) Bor a standard double diffused npn tran- 
sistor the net intpurity concentration is given by 


N(X) 


•“ Ng e 


-(xAg) 


+ N. 


-CxAg) 


B 


- N 


epi 


(2 39) 


where M(Xjg) - == 0 defines *ihe metallurgical jixnctions 

Unfortunately the use of (2 39) does not lend itself into a 
single analytical expression for our subsequent calculations 
Anyway by appTOximating the base in^urity profile by an 
exponentially graded profile 


N^(X) 




(40) 


a simple analytical solution is possible 
and (2 40) we match nCXj^^) * 


Comparing (2 39) 
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with 


^ jCX) ^jBO ^ given by 

n » ln(N(XjgQ)/^<XjQQ)) 


(2 41) 


and NCXj^q) and N(Xj^) has to be calculated from (2 39) VJe 
also note from (2 22b) 7} = Following a derivation some- 

what similar to that of Scott and Roulston [3B]v7e get better 
aporoximation of the ratio 


[^3 s 


<1 + a) 1 - 

T30 




. / OC, \ 
+ <7 ^ — ) 


«b®c Qbo 


1 + (1 + a) 


(2 42) 


where =« value of the ratio to be 

substituted by the left hand side of (2 33) or (2 34) as the 
case may be All other quantities are similarly defined 
Note that the ratio O^q/O^q so obtained is less than what wa 
have got previously without considering mobility variation 
The results thus support our explanation already giv n at the 
beginning of this section 
2 3 3 Variation of and with 

As already mentioned the OP model is based on the 
modelling of total majority carrier charge in the base region 
By assuming qua si -neutral condition excess majority carrier 
charge has been equated with injected minority carrier charge 
v^lle the latter is related with the corresponding current 
by charge control theory Thus the fourth and fifth terms of 
(2 9) (B^£l£ and respectively where and has the 

physical significance of forward and Inverse recombination llff: 
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time and B is a factor to Incaicporate base push out effect at 
high injection current) are arrived Guintnel and Poon accepted 
them as model parameters but as we shall see shortly they 
indeed show significant variations with collector current e g 
varies rcughly from w^/ to es the system passes 

from low to high injection 

In order to incorporate the effect of grading factor 
as well as variation of with we can model the excess 
charge in tejrnia of the excess charge density itself or by 
deriving an analytical expression of and as a function 
of corresponding current and still using charge control theory 
The first method has merits in a way that it allows us to 
incorporate the emitter crovxilng effect in a straightforward 
manner we will discuss here both the methods 

For its relevance we rewrite the expression (2 2a) 

here 

■ VJp " ^ 


and then use quasi-neutrality condition p ^ n + N^(x) 
Binstien relation kT/q and neglect majority carrier 

current j to obtain 


j 


n 


^^n^"irTS][j[xT aRc. ■** 


n 

n + 



(2 43) 


Ihis expression Is quite valid in normal operation even at 
high injection But in deep saturation and Inverse active 
high injection condition the negligence towards j^ would not 
be a fair deal since the collector region Gurnmel number is 
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far less than that of the emitter region and also because of 
capture ratio (1 e ratio of the capture cross section to the 
injected cross**sectlon ) is less than one for inverse injection 
the inverse active gain h^^,^ would be much less In spite of 
this fact since the current results from a short base diode 
and Jp from a long base diode remains far greater than 
Further owing to the low doping the epi collector region 
goes earlier into high injection than the relatively higher 
doped base region so we should expect h^^^ to Increase with 
^CE ™°'^^’^ou3ly until hi0i injection in the base region sets 
in But the situation is much mare complicated because of the 
presence of collector edge cixiwdlng (equivalent to emitter 
crowding effect for forward mode operation) the transport 
factor severely reduces which has an opposite effect of reducing 
hpgj ihe situation is thus complicated and to get exact 
analytical solution would be difficult one Nevertheless 
Gummel«>Poon model can be used though not rigoa^susly 


(A) Low Injection 

If everywhere n(x) «Nj^<x) equation (2 43) reduces to 

1 A 

■^n * ^^n N 3x ^*^A^ Multiplying throughout by Nj^{x)dx and 
A 

substituting the value of N^(x) from 40) the above e 3 q)re- 
ssion is integrated over the limit 0 to x Then a little 
rearrangement will lead one to [27] 


n<x) 


n<0) e 




qPft T) 


il5cAg^ 
{e ^ 


- 1 } 


(2 44) 


so with f » 



’7 


hn ^ ^ 

n(Wb) - £n(0) - ^ ^ - 1 ) 

Replacing by ( 2 44a ) gives us 

Jco = r^i [^(°> 


(2 44a} 


(2 45) 


On multiplying both sides by A_ the emitter area and using 

C2 

the normalised form of carrier densities (n ® n(o)/^.{0) 

O A 

« n(wj^)/N^<0)) this expression reduces to 


^cc “ ^fb<«o - ^ "b5 


(2 45a) 


with 

go £ 

^£o “ r^^a <2 46> 

As depends only on the junction voltage and only on 
^BE we will split the main current as so th t 

(2 45a) gives 


= Ij. 

f o o 


” I ^fo’^ 


(2 47) 


from (2 44) and (2 45) we get the following expression for 
excess carrier charge density 

t)x/wL n(wL> - fh( 0 ) i)x/v^ 
n(x) = n(0) e ^ [e ^ 1] (2 48) 


This expression can be utilised to get the expression for 

total injected minority carrier charge C^j|^ “ Q A^ / n(x)dx 

o 

Obviously this can also be split into a forward charge 
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depending on and thus on v^j, and a reverse charge depending 
on via Tiy. They are 




~ 9b< 


1,.(2L- l?f ,+ I 
(£ - 1)^ 


f(f -1 - 


(£ 1 ) 


(2 49) 


Here is the fixed base charge 

« Ae ^ <3x - q M Wb 

o 


f 1 X 

r\£ 


(2 50) 


(B) High Injection 

Now n(x> » N^(x) so (2 43) alTOplifies to 
5^ tt2qD^(dn/dx) This equation in conjunction with (2 40) 
gives the sinple solution for current density and carrier 
density 


•^n “ -^cc 


-S {n(0) - n(t^) ) 


<2 51) 


n(x) « n(0) * {n<0) - n<v^)l ~ 

In a similar way for low injection we arrive at 


(2 52) 


“ ^£o ^ ^ f ^ ’^o 


T - T 2 f - 1 n 

^r " ^fo t) ! % 


<2 53) 


Or “ %oll4T:"b 


(2 541] 
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Here again the current and charge are expressed as functions 
of normalised carrier concentrations To relate these densi- 
ties to Vgg and respectively we use the pn product at the 

junctions pn « (n + N^)n n^i^ exp(v/V^) where the neutrality 

condition is exploited In normalised form \ie get for n 

o 

and nj^ 


+ ”o> 










(2 55) 


Here 


ig is the familiar saturation current 
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f 


q D„n^ 


T "TC 
^ M, 




= 




°b< 


(2 56) 


Note that these set of equations are totally in agreement with 
Gumroel equation At low injection neglecting n^ in coit^- 
rison with 1 or l/f respectively we get from the sets (*5 55) 


a ^ the usual 


and (2 47) e ^ and 1^ = 1^ 

expressions of forward and inverse currents And for high 


injection n^ 1 then from the sets (2 55) and (2 53) 


^f == ^o 
^ ^ 2 
^o n 


V /2V 

■ 3 -fo ® ^ and = a 


ho *^s fo 


V2V_ 
, BC T 


(here 


(£ - l)/f) same as can be obtained from Guroroel 

equation 

Now comas the diversion from Gummel form Gummel 
has defined Knee current as the currant when minority 
carrier chairga equals majority carrier charge or 
from charge control theory bet us define as the current 
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where the two asymptotio ciarrents (2 47) and <2 53) intersects 
Thus 


vrn- 




or 


"k = Is « 




(2 57) 


V. Is the corar^sponding knee voltage Eliminating V. 
^ k 

from them 


“ ^olfo = (2 58) 

tJL 

Thus comparing (2 57) and (-> 58) we see = •s^ where 

2 

Tq ® wj^/2D^ is the simple forward transition time without 
grading effect At high current is reduced by two as the 
diffusion constant effectively doubles because of increase in 
conductivity in the base This result is also directly 
evident from (2 43) If we substitute n(x) -N^(x) there 
we would get = q(2Djj) ^ This is the san e^£pre 3 Sion as 

in low injection case except that is replaced by 2D 

From (2 47) (2 53) and (2 58) we note that at low 

injeeUon » ®ho^lk/®ho^"o injection 

“ ^fo^'o^ho * already defined or in 

general ^ ® ^^®ho ”o ^ 

g « 1 for n^ » l 

we now assume that g behaves monotonically with 
respect to the asymptotes of Figure 2 3 as a function of T) 

A fit function is proposed satisfies all these require 


ments 
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2 + a- ^ + 4n 
no o 

^ ° ^2 + a^Q)aj^ + 4n^ 


(2 59) 


So the general model expressions for currents and charges 


becomes 


^ ^ho hi 

‘2 + ^ho^^ho + ®ho ~ ° 


- ^ "fa T 

" ®ho •" "b ^ho ^ ' ' 


(2 60) 


2 + ■*• " 


(2 + a. „) ^ + n 
°ho o 

* % JL « 
>o x/f aj^ + "b 


... ®ho 


(2 61) 


tit I- 


^lo £ (T)^-"l)f + 1 ^tb * 

Iq Is subtracted in order to make If ® ^BB ^ 


ir - - lo for Vg^ « 0 

So far we have neglected the bias dependence of the 
fixed charge and the base width If we sredefine these 

quantities for zero bias they must be replaced by q^O^Q 

in the previous expressions \^en bias is applied where 
qj^ « 1 + A further consequence of these is that 

Ij^ and to be divided by the factor (see (2 58) and 

(2 56) respectively) These replacements incorporate base 
width modulation effect on the stored charges and currents 
and also on transit times An approximate eaq>ression of q^ 
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^ ""Tl 

would be (2 - e - e V(1 e ) where Cj^ is defined 

in (2 22b) For inverse mode operation Cj^ has to be replaced 
by Cp (» rj/w^) Also t) has to be replaced by x^/w^ (* Cp5^) 
wherever it arises 

Equations (2 61) and (2 62) are not real solutions but 
they give a realistic picture of the current and charges that 
obey (2 43) in the asynptotic cases Ihe advantage is that 
they enable us to derive a closed form expression for and 


X 

f 



X X 
£o o 


and x^ 
r 




% T 
ro o 


(2 62) 


2 

where we have taxen^^^ - Somewhat lengthy e3q>res8ions 

are available from (2 60) to (2 62) each for and which 
reduced to 1/q at high injectign and exhibit the limit 
(tj, - 1 + at low injection (idiere t) is positive for 

^fo negative for T^) 

As already mentioned §t the beginning of our 
derivation the variation of tJ^'anslt time can be Incorporated 
in two ways If we want to stay with GP model we can areuse 
(2 9) as 

Q 

Q, K a* cr^ + B i«. + T i (2 63 

^ ^Do ^1 fd f ro r 

wht^re we have used equation {2 58) and definition of and 
replaced and by their normalised form i^ (» 
ij. (® Further with help of (2 57a) equation (2 8) 

can be bought to the contact form 



(2 64) 


’GC 


- e /C^ 


where * ^cc'^^k normalised llnJcdna current 

Note that the model parameter used here through 
** ® corratant of the device structure It can be 

hypothetJcally defined as the knee current ^en there is no 
junction space charge effect whereas the actual Ij^ « 
is a variable quantity on account of Early effect 

The second way to attack the problem is to directly 
use (of course with prior modification) 0^ and from (2 61) 
Hence instead of (2 63) we should get 


2 + 


% 




+ n. 
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^ "b "b 


(2 + + n 

®ho ° 


q-a 


l'*ho 


I “ho ^ "b 


(2 65) 


\^ere and nj^ are to be evaluated from (2 55) 

So far at low injection we have talked about all the 
major effects excepting emitter crowding effect which we feel 
may not be so sl^if leant we will now consider high inject- 

ticjn effects then emitter crowding effect and conclude the 
chapter with a discussion on parasitica 
2 4 HICW INJECTION EFFECTS 
2 4 1 Conductivity t^odulation in the Base 

The onset of high inlection is said to occur at a 
value of current at which the minority carrier current density 
injected into the base becomes equal to the equilibri'um 
majority carrier density This takes place t^^^pically for a 
forward bias of about 0 7 v It should be realised that this 
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criterion does not represent a threshold beyond which there 

is an abrupt transition to high injection occurred Rather 

it should be looked upon as a limit beyond w - h the effect of 

high injection phenomena on the behaviour of the transistor is 

significantly manifested Actually the passage from low 

injection to high injection region is depicted by the smooth 

transiticoi of junction voltage dependence on transported 
V /V V /2V 

current from e to e in normal mode operation or 

e to e in inverse mode operation 

The high injection level effects render many of the 
assumptions and approximations in the low injection range 
invalid But the theory developed so far is quite general 
and can safely be extended to this region The only difference 
is that we need not require to use the rigorous expressions 
of and and also as the formers are insignificant 
with respect to and the last one quickly approaches its 
limiting value But some new effects like base push out 

emitter and collector crowding resistive drop in the base 
and collector region have to be incoaeporated 
2 4 2 Base Push Out Effect 

Let us try to make a clear picture of this effect 
before going for any derivation When the collector current 
is increased for a given collector base (^^ 3 ) reverse bias 
the concentration of carriers injected from the forward 
biassed emitter junction into the oollector base space charge 
region is increased Ohese extra charges add to the space 
charge on the base side v^ile they neutralise part of the 
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fixed inmobile space charge to the collector side of the CB 
roetallurgical junction If the minority carrier charge 
density is comparable to the doped charge d n ity then this 
causes the field distribution in the collector space charge 
layer to change That is the field on the base side increases 
While that on the collector side decreases It is the latter 
effect which has the dominant influence since most of the 
junction voltage appears across the collector side of the 
space charge epi-layer The fact that the total voltage drop 
across the collector-base space charge layer has to remain 
constant causes its collector edges to move further into the 
collector quasi neutral region to compensate for the reduction 
in field At the same time however the increase of f onward 
collector current results in an increase in the ohmic drop 
across the undepleted part of the epl-region caxising a reduc- 
tion of the available voltage bias across the p v junction 
Consequently the depleticjn layer tends to shrink one of 
these effectswill dominate in any given sit\iation However 
it can be shown that the end result is always an effective 
widening of the base width when the collector current density 
exceeds a specific critical value for a given reverse ^BC 
This is accompanied by an effective increase of and hence 
a fall in the common-emitter current gain 

Papers [7 14 16 ] dealing with base widening mecha- 

nisms generally make a distinction between the situation 
where the behaviour of the collector region which is next to 
the induced base region is ohmic (quasi sattxration) and where 
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It is space charge controlled (carrier swamped collector region) 
A unified theory has been developed by Rey et al [34 ] which 
works for any bias condition we will now prefer t a brief 
account of the work because of its relevance here and finally 
mould it according to our needs 

Figure 2 4a shows the situation when base widening is 
present The Influence of the electric field on the carrier 
velocity has been taken into account by the following pres 
cription 


U E 

when B ^ 

(2 66) 


when E^ < E Sg 

(2 67) 

>^0 = ''s 

^en B 

(2 68) 


where ** low field mobility Vg - scattering limited velocity 
are constant critical electric field For silicon 

7 

* 1400 cmA* sec v_ = 10 cro/sec Referring to Figure 2 4a 

O Q 

when the device is in quasi saturation state the collecting 
zone is an ohmic zone such that « qnv and 
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+ V, 


BIG 


(2 69) 

(2 70) 


vbere we have neglected the voltage drop in the quasi satu- 
ration zone (which is ro\ighly 100 200 mV) For E upto 

(v jC Br. n is equal to N_„ so that 

^ O op 




^ ^ep 


V 


(2 71) 
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VJith being kept constant If Is Increased we see from 

(2 71) and (2 66) that the carrier velocity v and hence the 

electric field E in the collecting zone groi/s This in turn 

leads to an increase of in accordance with (■? 70) It 

reaches E* for - q v_ Beyond this limit i e when 
s ^ sp s 

Jr* ^ ^ N v_ even if all the majority carrier moves at their 
c ep s 

highest velocity this would not be suffici^ant to carry 
It is therefore necessary for extra majority carriers to be 
injected into the collecting zone They are supplied by the 
emitter So above this limit one should use instead of (2 7l) 
the following 

Jc " 

In this condition a space charge 

P(X) - a (N^P - j^/q Vg) (2 73) 

appears in the collecting ssone this manifests itself in the 
modification of (2 69) to 

i = I tNep - ''s> 


Thus E(x) is an increasing linear function of x These two 
cases are shown in Figures 2 5(a) and (b) irespectively 
The description just given above assumes implicitly 


that Vqm j 

BE BIG s c (Figures 2 6(a) and 

Let us note the changes in the situation^|[by increasing 


while keeping j^ constant By virtue of (2 69) or (2 74) 
cG G 

as the case may be ds/dx does not vary and the area 
bounded by the curve E(x) versus x increases in a way dictated 
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by (2 70) Two oases must be differentiated (a) when ^ 

q N ^ V and when (b) j- > q v Depending on these 

(2 69) or (2 74) has to be used In both cas ^ ^ increases 

when Increases Here again we have assumed 

remairs Xorrer than E_w^ 

S ep 

rlnally the evolution of the distribution E(x) for 
Incr as^ng at constant shown In 

Figure 2 7 

ANAXjYTICAL solution so far we have seen different bias 
situation under v;hjch w occurs Let us derive an expre- 
ssion or it at different regions Using the appropriate set 
of aquitJons from the groi:p (2 66) to (2 74) the following 
relations are obtained for the current induced base width (w^ ) 


(a) 


(b) 


(c) 
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(2 75) 
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<1 «ep 1^0 ®C ^ Jc ^ •'ep ''s 
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WCXB = V,, 


t ^ “ep ''s 




-CIB “ “ep tl - ■-■aXT^ "" - ev^Ef 


C ‘ ep S 
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in 


(2 77) 


The critical current density designates the onset current for 
w'iich tie base widening phenomenon occurs and it is a function 
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of the 
(a) 


applied voltage (see Figure 2 8) 
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(b) Eg 
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(2 79) 
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(c5 V + V > E w 

"'CB BIC " S ^ 
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ep 
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•The base pu^ out effect Is Introduced in the expression 
(2 64) by introducing the factor B In the fourth term vftiere 

^ ''tiB ^ 3 

3 = ( -^ ■ % ) « (1 "“ T ') (since ***£- or In the fourth 

term of (2 65) by b « ** 1 - (since ^ 

we will replace by l 

A few points can be made In this connection Ihe 
actual picture of base push out Is much more complex than v^at 
Is depicted here Ihere are many side effects which render 
this derivation slightly inaccurate e g because of edge 
crowding the effective collector area over \AiiGh charge flows 
Increases and this leads to an effective increase of by the 
factor 


«b 


[1 


-4 


(i 

^o*E 


2 ^ 1/2 


(2 8i) 
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secondly when q ^ ^ approxi 

nation of Jq by (2 7X) is not justified to some extent 

Thirdly in quasi saturation state the oltage drop 
across the induced base region is though small but comparable 
to that across the ohmic region for a standard n^pin*^ structure 
ihis has to be taken into account Finally we have also 
neglected the zero bias depletion width when comparing 
and 

2 4 3 Quasi Saturation Region in the Collector 

L#*t us examine this region more critically In 
Figure 2 9 is shown the low versus characteristics of 
a typical high voltage transistor Here for a typical base 
current the characteristics exhibits saturation from 0 to 
A with a resistance equal to the n collector (Region X) 

From A to A (Region XI) the slope changes and the device 
enters its active region at A I'rom A to B the collector junc 
tion is reverse biased (Region III) At the point A th equa 
tlon is a ^ep^o^e region A to A is called the quasi 

saturation region In region I the p-v junction (we denote 
the epl n jregion by v ) is heavily forward biased which 
ensures the v region to be filled with injected boles with 
a concentration well in excess of the background in addition 
charge neutrality requires that the electron concentration in 
this region be appro jcimataly equal to the excess hole concen 
tration Thus the region is heavily conductivity modulated 
and supports a negligible potential drop 



Fig II I Output characteristics ( saturation) ot 
power transistor 



Fig 12(a) 6t Carrier concentrations in the P-»region 



Between A and A the P region is conducti ity modu- 
lated over a part of its length the remaining part will appear 
as ohmic resistance in order to model this ^'ejJon it is 
advisable not to use second equation of the set (2 55} since 
is not easily expressible rather one should equate the 

pn product on either side Defining = p(0)A.^ we then have 

o ep 

1 N,_ 

nj^(| + nj^) » (-^) Pq(Pq + 1) (2 82) 

In general is a complicated function of and l^ Let us 
derive a handy expression for p(0) The electron concentration 
in the v-reglon is given by 

n a p + N 

The electron and hole current densities are 


= q tip PE -qPpii (2 83)| 


The total collector current density * “^nc^ ^ 

negative quantity for the sign convension we have used The 
electric field in the collector region 


B(x) 



(2 65) 


has two components The first one is because of forward 
injection of the p-l> diode the second term is because of 
Ohmic drop across the epi -region Since our interest is in a 

regime where the transistor has a significant current gain 
then to a good approximation we should have j a* 0 and then 

pc 
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V , 

E(x) * ^ with dn/dx * dp/dx we see 


Jc ® -^ne * 2qD„(l+^)f| 


The solution of this equation gives 

P(x) » P(0) - ^ + -|® In ^ 


(2 86 ) 


If p(x) » the last term of the above expansion ean be 

omitted and p(x) Is thus linearly decreasing with x over that 
region This Is valid upto a point x » w^jq where 


'^CIB 


2q p(0 ) 


(2 87) 


(2 88 ) 


substituting this back In (2 86) we get 

A closed form solution of this equation Is not possible^ 
However one can approximate quite reasonably pCw^^j^) - N^p 
The potential drop over this conductivity mocul ted region Is 
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o ^ p(0) 




ep 
(2 89) 

The potential drop across the unmodulated part of the v region 
^ ^ ^ Ohmic and Is given by 

J/n ( W.„ ) 


o « /c^-epL.. 

q Hep 


(2 90 ) 
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(2 91 


vAiece 



/ 


The appropriate value of from <2 66) to (2 68) has to be 
substituted roughly the collector emitter voltage drop 

(Y^g) The termination of the quasi-saturation region for a 
given ^CE (or v^g) is given by the critical current (vide 
equation (2 79-)) 


The injected hole concentration p(0) can be determined 
from charge control considerations we have seen from (2 86) 
the dharge distribution is linear upto VJrlting as 

the hole life-time in the epi-region the stored charge is 
given by Q such that 


. ^ ^ qp(0)Wg^B 

■^pc ep 2 


(2 92) 


v^ere A is the cross section area of the epi-region Since 
charge density falls linearly with distance in the quasi -satu- 
ration region one can write 

i 




dx 




CIB 




(2 93) 


n 


so that from (2 92) and (2 93) 


p(o) 


Jc Joe ‘'eu 

2(.. V — |£ ~.gB) 


q 0 


n 


and 


D_ 1/2 

w *> 2( PP i 

”CIB J ' ‘ 


(2 94) 


(2 95) 


where modulation of the p-base region has been ignored here 
The onset of quasi saturation occurs at A writing 
this current density by j,j» and substituting x = w^ 
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2q D p(0 ) 

Jt = (2 96) 

ep 

The study of quasi saturation Is Important for a 
device acting as a switch since It Increases the turn on time 
and hence paver dissipation also 
2 5 EMITTER REGION SOLUTION [29 38] 

The doping of the emitter region is such that the 
semiconductor Is degenerate so far we have ignored this factor 
which Is O K for low level operations A closed form solution 
taking the variation of doping level emitter degeneracy band 
gap narrowing high level of injection etc into accotant is 
too much involved Roulston at al has derived the following 
rigorous expressions 


n(x^) « 






+ - Ng<X^) 




(e ^ - 1) 

•(Ng(Xj^) - Nj5^(3<^))e + Ng(x^) *" 


(2 97) 


(e ^ - 1) 


(2 98) 


where n(3^) p(x^) denotes the minority carrier charge den 

sitles at the base and emitter side of the junction depletion 


region potential across the junction and 

V — * V V 

^BIE ^BE 
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2 6 EDGE CROWDING EFFECT 

The analysis of the transistor is mainly based on the 
one dimensional current flow assumption Ho /ever In any 
practical transistor there will be a cross flov; of base current 
to ;ards the base terminal This flow of base current parallel 
to the emitter base junction will cause a transverse voltage 
drop which makes less of the applied base emitter voltage to 
appear at the centre of the junction as against the edge (see 
Figure 2 11) This increases injection of minority carriers 
from the edge of the emitte'^ as compared to that at the centre 
thus causing the current to crowd around the periphery of the 
emitter (close to the base metal contact) Emitter edge 
crowding becomes significant when the difference in voltage 
between the centre and edge of the base is of the order of 
thermal voltage While studying low Injection operation 
the above phenomenon may be ignored as the base current Itself 
(and hence the lateral Voltage drop) is small However the 
currant crowding at the base region underneath the edge of the 
emitter causes earlier onset of the conductivity modulation 
and other high injection j^enomena there The most serious 
concern with this effect is that it reduces the power rating 
of the transistor For compact modelling it must be taken into 
account though that leads to unavoidable complicacy 

Several papers [12 16 23 25 ] are listed at the 

end which deal with this effect Owing to the inherent con^~ 
lexity of this problem mixed with other effects everyone 
has to make some ridiculous ( ) assumptions to get closed form 



Fig 211 


Fig 212a 


Fig 212b 



Fig 2 11 Typical cross section of a transistor witf 
transverse axis(l)chosen 

Fig 212(a) Showing the voltage drop occurs acr^ 
the base under different approximatl 

(b)Comparison between actual voltage dJ 
& that nrtCilMi<-*«/»-4 -/oAft ^ 
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analytical solutions we strongly feel that these does not 
give true picture of the situation over most part o£ the opera- 
ting region As for example the papery [23] by J R Hauser 
seems to be most reasonable amongst them But it starts with 
the assumptions amongst which (a) the ignorance of conductivity 

modulation in the base v;hereby R. is made constant and 1 

v(x)/v ^ °° 

goes as i^ e (low injection limit) (b) is a constant 

s 

are most seriously objectionable especially in high injection 
regime where this very effect is important 

we propose to give some new model under most reasonable 
assumptions They are 

1 The emitter region represents an equlpotential plane 

2 The vertical plane perpendicular to 1 axis at 1 « 0 the 
edge of the emitter contact is at equlpotential 

3 ^ ® function of voltage (this may be omitted) 

4 We start with single base stripe (I type) metal contact 
though the theory can easily be extended for double stripe 
contacts 

Let us start with seeking how the potential profile 
changes as one passes from edge to the centre in the base 
region Throughout our discussions we shall assume the resis- 
tance Rjjj from the base metal contact to the section (1 « 0) 
in the base is a constant quantity Xf we can replace the 
part of the base region beneath the emitter region by a lumped 
resist nee and further assume the whole of the base hole 
current traverses this distance (Lg) being originated from 
the centre of the Junction (point injection) then voltage 
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would have fallen linearly (for rvgti transistor) from edge to 
the centre with a gradient so that V(l) - 

-rjjljjdAig) (see dashed line - - in Figure I2a) Because 
of distributive effect the gradient (magnitude) itself rather 
linearly decreasing and finally becomes zero at 1 = Lg We 
then vxjuld have V(l) = - (s*5e dotted 

line in Figure 2 12a) where we have assumed the base 

current density Jjjd) is uniform) all over the base region 
But this is not the case as the voltage at the edge (1 - 0) is 
higher than that at the centre (1 * Lg) and cvirrent density 
has exponential dependence on the base emitter voltage (the 
emitter surface is equipotential region) Hence larger portion 
of is originated from the edges th n the centre this fact 
brings about a third order term ( (1A>) ) in the expansion of 
v(l) (continuous line —— in Figure 2 12a) and so on The 
conductivity modulation in the base complicates the situation 
in that case the coefficients of (l/h)^ in the expansion of 
Vd> suffer change For high base currents the base profile 
shows sharper fall upto L followed by a slow fall which ends 
at Lg (shown by continuous line in Figure 2 12b) In that 
case expanding v(l) in the pcswer series of 1/Ijg would not be 
justified we propose to use a scheme whereby V(l) has to be 
modelled by a series esgpansion of (X/1 j) upto h (dotted line 
in Figure 2 12b) and then by a constant value (dashed line in 
Figure 2 12b) from L to Lg The point L would be determined 
by the fact that the positive and negative error (shown by 
hatched region) cornea equal 
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We now Introduce the following dimensionless variables 
X * l/li normalised voltage drop v<x) * 

v(l) * V(l)/v^ in the base According to th^» aoove scheme 
we postulate 


v(x) 


CO 

S 

n=l 


C ^ 


for 0 < X < 1 


S v(l) for 1 < X <1/1^ 


(2 99) 


where s are parameters (function of 1^) to be determined 
Mormally upto thijcd or fourth term is retained in the expan 
sion of (2 99) which ensures sufficient accuracy F»r latter 
convenience we put down here the following expansions 


^v(x) 


B (n + for 0 i x < 1 

n=o 


v('<)/)l|g 

e 


2 (n + l)k_x 
n-o 


(2 100 ) 


where is another set of parameters Prom (2 99) and (2 100) 
one can establish the following identities 

i H ’‘a ” 

(2 101 ) 

» •|(C 4 + C^/2 + Cj^Cg + and so on 

ihe parameters k^^^ can be obtained likewise by replacing 
by C /n in the identities of (2 101) we further use the 
following notations the sheet resistivity in the base 
region defined as the resistance of the slab of cross sectional 
area (wiare Hg is the breadth of the emitter junction) 
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and length unity Using (2 40) and (2 50) we see its value 
given by 


‘ba * 


(2 102 ) 


where is the same as v/hen no bias is applied Hence 


TkCx) 


bo 


l/(ttp5^(x)/Lg) = ^ 


(2 lO^a) 


(In absence of conductivity modulation <2 

At the emitter edge in the base (x - 0) we have 
rj^(O) = From intuition we can at once get the value of 


Cl 3S 


= dv(x) 


35 — U ■ 


k-ffliil . 1 


r 


(2 101 ) 


Let us divide the junction surface area into a large 
number of stripes of length Hg and thickness dx The base 
current density jj,(l) would be uniform over any typical stripe 
and Varies monotonoxisly from stripe to stripe an view of 
(2 99) The current originated from a typical one would be 

K 

jjjCl) Hg dl and the total current ^ ^ 

1^0 

Transforming the integration variable from 1 to x along with 
the help of (2 13) we can write 


« 1 . 




''be<o> 1 


/ dx + 


e 




y 

Q 


1 v{x)/f) 


dx 


(2 104 


where * Vgg(o)/V^ This on inserting the expression 

for from (2 100) transforms to the following identity 
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1 = 1_ [i, >t_ + 1- E k ] 

r *• lo n 2o en 


n=o 


n«o 


where 


t a; ^ 0 

to Ib ® 


"20 - i; « 


(2 104a) 

(2 104b) 


One can take the integration upper limit in (2 104) to 1/lj- also 
(A) Low injection cage 

The voltage drop at a distance 1 from the edge in the 

base 

1 1 

V(l) = - / dtbd) [lb / 3b^^ ) Hjs dl ] 

a o 

with drjjj(l) = (ignoring conductivity nodulatlon ) and 

transfo arming the variable index from 1 to x 


v(x) =--^ J “ Vr ^ ^ 

» - 1 / dx [l - 1 (1, S k ix” 

r ^ *• r' lo n-1 

■^"2o„lj!^en-l=^l 

* ^1* ” ‘^l"r t"lo “ "20 ’*an-2^ 

(2 105) 

Note that we have debarred fjDom giving upper limit to the summ- 
ations It is an option to the user to specify it and normally 
4 or 5 is sufficient Comparing (2 99) with (2 105) we see 
for n « 2 to <» 

S “ ~ Vl^^lo’^n-a ^2o^en-2^'^” 


<2 106) 
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where C. Is given by (2 103) 

for 

The expression^co Hector current (density In normalised 


form (from (2 63) and (2 64)) /ould be 
Vg-(X)AT 

e -v(x) 


l^(x) ^ 


%orr 


- igCo) e 




(2 107) 


Now %(x) - ■*' ^2^^^ 

« Defining shown 

(2 107) to be approximated by 


ij,(x) e [1 +q,,.{l - e''‘-^h + 

(2 108) 

The corresponding expression for the Unking current (taking 
upto second term in (2 108)) 


n + 1 


*cc 


Vc<«> IJ, Vi * -K J. Vi« - “ im;* > 


* ^ (VI - 


%L 


(2 109) 


If we take into account the modulation of base reals*- 
tlvity then “ Cj^CO) [%(3c)/<^(0)]''^ where the bracks tted 

quantity is already evaluated while passing from (2 107) to 
(2 108) Equation (2 108) 1® niodifled to 


v(K) * 


- ^2r 


2 Va*" - <^lVlo t L 


n 


n=2 


n=2 


^ n»2 n % ^ 


n 


(2 110 ^ 
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Where we have neglected the recorablnation-*gene ration current 
component of as it is Insignificant in this and high injec- 
tion region Correspondingly (2 106) modifies to (for n L 2) 


'n 


C-lJl, k ^ 

1 r lo n-2 *■ n 


TJ 


^2r 


^2r 


n-i 


r lo 


S (n + l)k k ] 


(2 111 ) 


and would be given by (2 103) The sequence for evaluation 

of these parameters are as follows m take 1^ * 1 /1^<0) 

r cc c 

and then find out the parameters (with the help of (2 101) 

<2 103) and (2 111)) according to the chain ^x: ^ ^ "* 

C 2 ** ^2 upto the number that depends on the user ihen 

use these k^ s in (2 109) and solve it by Newton-Ralphson 
technique to get 1^ The parameters finally can be redeter- 
mined by using this Ij. only one iteration is sufficient as 
the final result dcies not depend critically on 

The dc or low frequency base resistance can be calcu- 
lated from considering the power dissipated by the base current 
Hence we define effective base resistance by 


% 



)Hg di f 


or 


1 0. 

r *^bo 

o 




K i, 

lo 


S k 
n 1 


h-1^ -* 


(2 112 

At low injection with no conductivity modulation the retlo 
rj^/rjj(x) can be taken as unity in that case (2 112) is 
itforked out to be 
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"b - Vboti - 2 w „!» j, ri 

n-z n 2 n «! 

)c k 

— n"TV^^ 

where designates the said resistance for low Injection 

line 

no conductivity ittodulation If one would like to incorporate 
the latter effect then appropriate ejcpression for the ratio 
^ Inserted as have already been done in 
(2 108) With this modification (2 112) is again solved for 
we here put down only the final result 


lioon 


"’line ^ 


SS^ 


n n — 1 n 


n 1 
n 


" ^r^lo 


n-X n-n 

Z Z S 
n«2 n n «l 


k .k k 

n-n -n -1 n n 


n 1 m 
n J 


(2 114) 

where is the base resistance and low injection with 

Ucon 

conductivity modulation Physical interpretation can be 
given for the occurrence of the various terms Jn ('> 113) and 
(2 114) We see for extreme low injection whence 1^ ® 1 
** 1 and all other ** 0 both the expressions reduc to 

rij/3 

<B) High injection case 

The approximation that holds in this region is 
^1 ^ perform the same series of calcu- 

lations as has been done in low injection case we will 
put down only the salient steps here 



59 


The normalised collector current density i (x) is 

c 

given by (2 107) We reformulate it as 

1(0) q. q, 1 


fo *^0 


or l^(x) e i e''M /2 + _2i eV(x)/ 2^~^/2 

° *^20 

K . /c‘°5'3bo 1/2 

where and = b 1 


(2 115) 


(2 115a) 


7 ^7 

Defining eqn (2 115) in expanded form 

^2o 

(upto second urder) appears as 


i^(x) B 1^ [a''('‘)/2 - q^ + ^ 

(2 116) 

with 


= l-o<l ” ^ 2 ‘Ih^ 

(2 116a) 


The base resistance as a function of x is given by 


rb(x) 


^bo 


(g-v(x)/2 . 




- wv(x) 
e ^ ) 


(2 117) 


Using (2 115a) and (2 116a) in (2 117) at x = 0 


we gat 


rb(0) 






(2 117a) 
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The expression for v(x) is derived with the help of (2 117) 
and the resulting expression is compared with (2 99) This 
at once gives an expression of similar to that (2 106) in 
low injection case 

n-l \ n-l ~ '^lo^r 


n-2 k Jc. 
{ 2 - n-n ■»2 : 

n-2 " 


' % 


n + 1 


n 




(2 118) 


where c. is as in (.'> 103) The coefficients k_ „ (where ») 

can be 1/2 -1 etc ) has to be evaluated from (2 101) in the 
same way excepting that s has to be replaced by 

The corresponding expression for collector current 

density 

1/lj. 

^cc ** ^r 

o 

on using (2 116) we finally get 

^eo “ ’^-1/2 n-l “ ‘^h ■*■ f *^1/2 n-i 

+ (^ - l)e''^^^/^ (1 + 2 ^ a v(l)/2j ^2 ^^g) 


The effective resistance at high injection 


‘hh “ [’‘-1/2 n-2 ‘^‘-l n-2 - Vlo 


("'s IL-t-i (k_^/2 ^ - <^k.2 „ )k 


n =o 


h-n 2 


(2 120 } 
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rb(0)l 

where extreme high injection 

^bti ® *'bh<°> Cl “ % ■ 2 Vio<l " ^2 120a) 

Thvis at extreme hi^ injection the base resistance £alld 
Inversely with the collector current 

So far we have developed the sit\)ation for single 
base stripe contact For double stripe contact the theory 
can be easily extended In this case l*g « half of the emitter 
length appc.ared in the previous expressions should be 

considered as half of the terminal currents At high Ihjec 
tlon the theory developed remains essentially the same At 
low injection one should include a factor (1 - Ij-X) to 
dr^Cl) or dr^Cx) wherever it arises 

The sheet charge density of excess carrier elso 

appears to be modulated by the edge crowding effect The 
expression for qb£ and is given in 65) and (2 55) respec- 
tively where is a function of x because of gw crowllng 
effect Defining 

^o n^(0)(t + n^(0)) (2 121) 

where the first suffix with n^(0) represents it is at the 
emitter surface and the zero in the bracket defines the quan- 
tity at X = 0 or 1 » 0 we at once get from (2 55) and (2 121 ) 

<2 i21a) 

00 o 


Then 
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o 

1 

= Ijf [ / %£(3<) <3bc + (1/ly - 
o 

This integration can best be evaluated by numerical techniques 
using Legendre-Causs quadratvtre rule The integration limit 
is chang'jd from (Oto l) to (-I to +l)by introducing the 
variable I so that x « ^/2 + 1/2 First n^(x) has to be 
evaluated from (2 121a > for a particular x and then to be 
used in (2 65 ) to find out q 5 £(x) This is then Inserted in 
(2 122) and the inters tlon is carried out as summation 
2 7 PARASITIC EFFECTS 

Schematic cross section of a typical transistor and 
its parasitic effects modelled by equivalent lumped elements 
dre depicted in Figure 2 13(a) and (b) respectively At 
high frequency the equivalent circuit is much more complicated 
as the reactance because of lead inductances also start 
becoming significant The diode l >2 there to represent the 
lateral injection from emitter to base through the side walls 
of the «nitter This is a long base diode Similarly 
models the diode action through the inactive part of the 
collector region which is Important when the latter itself 
is acting as emitter in inverse active mode This is a^eo a 
long base diode The capacitance is the capacitance 
between the emitter and base roetal contacts and 
oapecitance of the inactive part of the base region One 


can assxme 



^xiae 



(a) 
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(2 123) 


where r. Is the ratio of the active to inactive part of the 
base region Roughly 


1 - 


and the voltage is given by 


^be ’*■ " “^0^0 


(2 123a) 


(2 124) 


*1116 modelling automatically in ^be3 

In inverse active mode GP model miserably fails 
This is because the basic assumptions behind the derivation 
of model are violated GP model remains valid when a is 
close to 1 But in inverse active mode a differs much from 
unity Since the coll*>ctor doping level is much less than 
that of base or the ratio G^/G^ W is not large The 

collector (acting as emitter) injection efficiency is small 


Secondly because collector to knitter area ratio is less 


than 1 the transport factor is also considerably less than 1 


In oxrder to stick with GP model we once again start 
with tha basic aquation (2 2a) and divide throughout by 


pjip 5^ to gat 


n 

j„(i - 

h p lip jp 


d 


tin 


(2 125) 


and apply it over the intrinsic sytmatrical transistor 
(shown by dotted extension In Figure 2 13a) We cannot 
neglect the hole current J this time and write jp/3„ “ 

(%/Gc “ using quasi neutrality condition p = h + N^ 
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we reform the above equation to 

n ^ 

n + ~ ‘ p ' cE 125a) 

where injection n « at any point in 

the baee Integrating the above expression from collector to 
emitter assuming transport factor to be unity t/e arrive at 

Jn< j" ^ " ‘‘rSfae ^ f"P3? <2 lisb) 

3C^ n ap^ "n ^ 

Multiplying by A^qall throughout and as suiting to be 

constant 

^n^^ ^r^bc9br^ * 125c) 


we finally get 


2 n A 

. _ °n *e " ® 

” OboUi+Tl lir«oo^%r^ 

(2 126) 

lg(e ^ - 1) 

+ (l - Hr9bc>%r 

(2 126a) 

where is given by (2 8) and 

9 


, ; 

«bNA 1 - 

9bc “ “5 ST 

0. coth -S2 

®P P P 

(2 127) 


This is the expression for transported current in inverse 
active operation since the collector region enters rather 
earlier in the hi^ injection region than in the base this 
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nD„/(wj,N^) 
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(1 - 


£“^)( 


coth -|2 + r^Jg/q) 
ep p P 


(2 132) 


We thus see the gain in inverse active is much 

smaller because of low doping to the ooHsctor site (N^p << 

N^) and because of the presence o£ the overlap diode 0^^ 

(manifested via J^) 

s 

2 8 SUMMARY 

So far we take the pain to study in great detail 
about all the knovn major events that occurs 
in various modes of operation of a transistor It is not 
viable to incorporate all of them in a single package and 
extract the relevant parameters with acceptable accuracy 
Anyway one can cleverly pick up data from certain regions 
v» ere some or the other effects are predominant and hence 
the relevant parameters may be extracted which can be subse-* 
quently used for modelling other regions One may have 
to reiterate among the various stages until an universal 
convergence is reached In our discussion we have indicated 
different approaches for modelling a particular effect 
However we will choose the one vdiidh suits well within op 
model In Chapter IV we will first consider some of the 
nunnerical techniques and then we will make a systematic 
approach to extract the model parameters 
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CHAPTER 3 

TRANSISTOR CURVE TRACER 

3 \ INTRODUCTION 

*The modelling of higher order effects (eg variation 

of base resistance with collector current) requires accurate 

measurements and special set-ups Accordingly a special 

purpose curve tracer has been designed it operates at frequ 

ancles around l kHz so that both transistor noise as well as 

power line Interference can be averted The system is 

insensitive to the scanning input waveform Moreover it can 

be tBsed for measurements on JPETS t OSPBTS UJT etc in 

addition to BJTs 

3 2 BASIC BUILDING BLOCK 

The basic concept beMnd the system is formulated 

in the block diagram as shown in Figure 3 1 The basic idea 

is simple The input ac signal is converted to an ideal 

rectified output excepting for a possible multiplication 

factor It is then applied to the collector point via a 

small san^ling resistance and collector series resistance 

R if any The same input signal is sensed by a ZCD (zero 
c 

crossing detector) or a peak detector (both positive and 
negative) which gives digital output levels This is sub- 
sequently fed to a frequency multiplier (so that for any 
collector sweep a new base step is generated) The output 

signal is used as clock of ah up/down counter The digital 

outputs of the counter are fed to a DAC followed by t 
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Fig 3 la Block diagram o£ the curve tracer 

multiplying DAC (rather acting as an attenuator) so as to 
generate saw tooth/triangular wave as the case may be 
This output is then processed successively through a buffer 
and then a linear VCCS The output of the buffer (V) as 
voltage step or the output of the VOCS as current step can be 
applied to the input terminal (T' It Is then connected to 
the base (B) or emitter (E) while the otht-r to tlie ground by 
losing the ganged switch depending on whether wa are interested 
in CB configuration or CB configuration if the signals from 
appropriate points are picked up (© g in order to plot 
output characteristics i e 1^ family of curves with 

Ijj as parameter the c point is connected to X-channel and 
Y point to Y channel while keeping oscilloscope in XY mode) 
a steady picture would be available This is then stored for 
subsequent use or for getting a hard cojy 
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3 3 1 Collector Sweep circuit 

The system indicated in Figure 3 3 gives full vave 
rectific tion withA/ithout inversion and jith a gain 
controllable by one resistor followed by a current booster 
clrciiltry This has the advantage over ordinary rectifier 
that it does not suffer from zero crossing distortion If a 
sinusoidal voltage whose peak value is less than the cut-in 
voltage V (06V) of a diode is applied to an ordinary 
rectifier circuit we get output to be zero for all the times 
And for a low level signal there is appreciable phase lag in 
zero crossing By olacing the diodes in the feedback loop of 
an OPAMP the cut in voltage is divided by the open loop 
gain Ay of the amplifier - Vy/lAyl 60 ^iV for Ay * 10^) 

and the diode acts like an ideal rectifier Though a straight- 
forward modlfloation of the ordinary rectifier (whereby each 
diode is replaced by an ideal diode so described and a 
summer is introduced at the end) is feasible but the choice 
of our circuit is obvious as in the former case each input 
OPAMP is strained by large differential voltage across its 
inputs and also it requires an extra stage a summer at the 
output Moreover the gain is easily controllable in our 
adopted system by only varying the resistance R^^ 

The basic operation of the circuit is as follows 
Let us consider first the half cycle where is positive 
with the ganged switch position as shown (used for measurement 
on npn transistor) D 2 is ON and is OPP since cond 
ucts a virtual ground exists at the in|mt of OPl Because 



b) [Y;y[vVYVVYVYYYVYVYVVVYYVVYYyVYYYYY^ 





Fig 3 2(q) Input voltage waveform (b) Collector voltage 
derived from(a) Using ZCD (c) The same as 
(b) but using a 90 phase shifter preceding 
ZCD(d) Ic VcE characteristics with collector 
voltage os in(c) 


71 


3 3 1 Collector Sweep circuit 

The system indicated in Figure 3 3 gives full vave 
rectific tion withA/ithout inversion and jith a gain 
controllable by one resistor followed by a current booster 
clrciiltry This has the advantage over ordinary rectifier 
that it does not suffer from zero crossing distortion If a 
sinusoidal voltage whose peak value is less than the cut-in 
voltage V (06V) of a diode is applied to an ordinary 
rectifier circuit we get output to be zero for all the times 
And for a low level signal there is appreciable phase lag in 
zero crossing By olacing the diodes in the feedback loop of 
an OPAMP the cut in voltage is divided by the open loop 
gain Ay of the amplifier - Vy/lAyl 60 ^iV for Ay * 10^) 

and the diode acts like an ideal rectifier Though a straight- 
forward modlfloation of the ordinary rectifier (whereby each 
diode is replaced by an ideal diode so described and a 
summer is introduced at the end) is feasible but the choice 
of our circuit is obvious as in the former case each input 
OPAMP is strained by large differential voltage across its 
inputs and also it requires an extra stage a summer at the 
output Moreover the gain is easily controllable in our 
adopted system by only varying the resistance R^^ 

The basic operation of the circuit is as follows 
Let us consider first the half cycle where is positive 
with the ganged switch position as shown (used for measurement 
on npn transistor) D 2 is ON and is OPP since cond 
ucts a virtual ground exists at the in|mt of OPl Because 
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FULL WAVE RECTIFIER CURRENT BOOSTER 



To meosureme 
section 


To collector 


point 


Fig 33 Collector sweep circuit 



Fig 34 Differential voltage generator 
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is nonconducting ther^ is no current In the R which is 
connected to the non -inverting input of 0P2 Hence this 
input is grounded The system then consist of tvx3 OPAMP in 
cascade with the gain of OPl equals (- that of 

0P2 equals (-1) The net result is 




for Vjj|> 0 


(3 1) 


Consider now the negative half cycle of v 


IN 


This tinve D, 


is OFF while is ON The inputs of OP2 stay at the same 
potential v since the negative input of OPl is at the ground 
(virtual) potential the Kirchhoff s ctirrent law at this node 
gives 


IN ^ V ^ V 
■*' 2R R 


or 


" -1^1 ''in 


(3 2) 


And that at the node (negative input of 0P2} gives 


V - V 
o 


V 

15 


or 


''o' I'' 


^ ''iN^ 0 


for 0 


(3 3) 
we thus 


where use is made of the previous expression (3 2) 
see the outputs for the two half cycles are identical pro- 
vided the four resistances R are perfectly matched In actual 
practice using a variable resistance in place of R connected 
between snd negative input of OPl da enough to avert the 
problem that may arise because of mismatch of the resistors 
The effect of input dc offset voltages and currents is more 
or less nullified by using the method discussed in Appendix i 
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since the output voltage is always positive this can he used 
as col lector sweep voltage for npn transistor VJith the 
ganged switch pushed to the other position A and S 2 ** B) 

enables one to use the same system for pnp transistor as here 
one gets Inverted rectified output 

An OP AMP (we have used 7410 can provide output 
current typically some 20 30 rtiA as a roaxiinum The high injec- 
tion region of the transistor well exceeds this limit So 
what we need is a current booster after this stage The 
circuit is essentially an emitter follower with the opamp used 
at the input to get rid of the base emitter voltage drop 
suffered by the output voltage to be taken from the emitter 
point *1110 rating of this transistor should be greater than 
the rating of the transistor on which measurement has to be 
done one can use field effect transistor (FBT) as wall 
The switch (S) position is shown in Figure 3 3 while measure- 
ment is taken on npn transistor For its complementary case- 
study the switch position connects negative supply to 

the pnp transistor so that the BB Junction is always forward 
biassed and negative feedback through the OPAMP is fulfilled 
When the output characteristic (I^ vs plot 

keeping as parameter) is wanted the switc^i is to be 
connected with point B (as shown in figure) and then the 
collector series resistance Includes the sampling resistance 

R-j also For measurements like vs keeping C(XiS- 
s 

tant comes from powerr simply (negative or positive supply 

iLl V 

d«)endlnB <»» whethar npn or pnp translator la chosen) 
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the switch is connected to F H is sno’-te^ so that a 

c 

constant voltage appears at the collector oolnt whose mag'll 
tude can be specified by vaiying Since now comes in 

the feedback loop its magnitude must be taken mall When 
ever the next base steo ( Ai^) comes the voltag at the out 
put of the opamp or base of the pover traisistot uf e "S 
in bantaneous change of amount v ln(l R_p/l_ wher 

* o B 

pis the small signal gain Slew rate of the OPA IP ould then 
restrict this instant neous change and or la»-ge collector 
Gurr nc the out ut of OPA'IP may gat saturated unle s a mall 
Rg is used Tn best remedy -For this oroblem is to u»e a 
power FBT in tead of tne BJT and to pla^-e Rg net een tie drain 
point and the supply point since there is oractically no 
gate current the same oun^ent flows at tne dra^-n and sou»*ce 
♦'erminals and s again the (negativ ) input or the OPAIIP 
draws no app»*eciable cun.ent this is the same as the collec 
toi cui^'^ent itself 
3 3 2 Measu»"ement section 

The collector cu^ent is measured by samoling the 



itsel*F This is notiiirg but the emitter cuYent of the oower 
transistor as the negat ve inout terminal i not drawing ary 
significant cunrent o^nts are keot in lind /hile 

designing the in t'^umental amoli ter Fiistly it should 
not draw any appreciaol cu'-rent from eithe input and 
secondly the common mode voltage rejection should foe as far 
as possible The fir t ooint is obvious second point arises 
because of the fact that there is a large exertion of cominon 
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mode Voltage at the two ends of Rg whejreas difference signal 

happens to be very small at lower base steps One can use 

larger value of R- but that would lead to oth r complice 

s 

tions as has been just discussed Popular Instrumentation 
amplifiers take care of the first issue but falls to meet 
the second requirement 

A modified form is adopted here AS depicted in 
Figure 3 4 the input pair of LMllO (dedicated voltage foil 
ower) serves the purpose of a buffer stage Voltage offset 
balancing is done as shown One of the outputs is inverted 
so that the current due to difference signal alone is drained 
throtigh The output voltage is then 

V - RfA (Vj - Vj) = Rj/R Rgljj <3 ^ 


Hie output voltage Is thus proportional to with the 

choice Rf = R =• 10 00 K Rg = 1 000 K the magnituie of 
(in volts) gives the collector current in ttiA 


In choosing the OPAMPS LM108 is preferred because of 
its availability and its muicSh lower offset drift than ordinary 
general purpose 741 type of OPAMPS The necessary compensa- 
tion is made The circuit is tasted with a cjornmon mode 
voltage as high as 10 V but with an output of 2-3 mV only 
whereas popular design ensures the common mode output not 
less than some 30 mV or so under the same condition 
N B It is to be noted that since the voltage V^ Is itself a 
regulated voltage source the upper follower (bMllO) is 
redundant and it can be omitted as well In that case if 
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IjM 108/LM110 have good delay matching -V 2 and will be in 
good antiphase and hence this ensures better accuaracy 
3 4 BASE STEP GBMERATOR 

3 4 1 Synchronizer and Frequency/ Multiplier 

As has already been discussed in section 3 2a base 

step change should occur at the beginning of every collector 

sweep (1 e for = v ~ O) so that a nt,at display of the 

family of curves will be obtained But if we are interest««d 

in the saturation (low region than the noisy scene close 

ce 

to V ^ 0 that arises because of Imperfect synchronisation 

ce 

may not enable one to extract good reading from that region 
It is then advisable to make the occurrence of base steps 
v^en applied voltage and hence is at its peak The 
type of display we will 9 ®^ is shown in Figure 3 2(e) Due 
to charging and discharging of parssltic capacitances some 
uncharacteristic traces or loop would be formed at the fag 
end of the curves a region vAiich we are not interested for 
low 1^ measurement Anyway the modification that we 

require for the latter case is to use a 90 phase shifter 
proceeding the zero crossing detector ihis part of the 

circuit is shown in Figure 3 5 

With the ganged switch position shown (s^ 4 and 

Sc **2) the input OPAMP (OP4) is activated as an integrator 

with a transfer function given by [3] 

g 1 t 

V<s) - - Avo 

Vijera Is a negativa number being tha open loop voltage 
^ ain and Si tlw dominant pola of the OPAHP in absanca of c 
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The transfer function has two poles on the negative real axis 
as compared to a single pole at the origin tor the ideal 
integrator we observe that the performance o the real 
integrator departs both at low and high frecjuencles At high 
frequency the integration is limited by the finite bandwidth 
(- Sj/2 7t) of the op rational amplifier while at low freq 
uencies the integration is limited by the finite gain of the 
OPAMP For pA741 A^^ » 10^ f^^ » 10 Ha and since our f 
is limited between 200 Hz to l 5 kHz the expression (3 5) 
reduces to 

hence an ideal integrator The choice of RC is defined by 
the attenuation wanted With the same resistance R applied 
to positive input only input offset current flows through 
capacitance The resistance R itself has to be chosen low 
enough so as to mask the effect of offset current and offset 
voltage we have cjhosen R « 10 k C =» 0 05 RF so that 
Ia^^ I =8 2 for f » ^ 636 Hz Its variation between 1 to 3 

(roughly) for the operating frequency range 300 to 1 kHz 
IsaUright as we are interested only in the phase shift and 
sitosequent stage have nothing to do with the magnitude of 
the voltage 

With the ganged switch rotated to the other position 
(S^**l 3 g'» 2 ) the integrator circuit is isolated and it 
is in fact modified to a mere voltage follower output always 
sitting at zero volt because of feedback action through R 
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Now Is directly connected to the next stage a zero 
crossing detector Since we aim to get pulses at both posi 
tive and negative going zero crossings e <^00 s LM361 It 
is a high speed differential input comparator [ 1 ] with comp 
lementary TTL output voltage levels It provides Independ nt 
strobing facility for thi- complementary outputs and can be 
operated over a wide supply voltage range It giaaranteas 
low input offset voltage and tight delay matching on both 
outputs 

The complementary outputs are Individually sent 
through a simple RC differentiator and then thr ugh a diode 
to clip off the negative spikes The value of R-is chosen 
so that the spikes are sufficiently high and the selection of 
C is dictated by the demand that the time constant RC should 
be small so that the tail drops close to zero before the 
next spike comes The choice R *» 51 k and C ® 1 pP gives 

good performance The outputs of the diodes are then sent 

through Inverters to piirl^^ the pulse shape and then through 
NOR gate Pulses are interdlgited at the final output the 
latter thus has frequency twice the frequency of the indl 
vidua 1 ones 

Por the two inverters and the final NOR gate three 
fourth of a 74C02 (quad two input CMOS NOR chip) is used 
Here choice of CMOS gate is essential as with input reals 
tance 51 k the input of TTL inverter gate would never go 
to lower logic level improper choice of R and c in the 
differentiator circuit aging or hardware troubles may 
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distort the shape/levels of the output pulses To rectify 
this a level detector can be inserted with the other input 
of the comparator maintained at 2 4 V 
3 4 2 Trlangular/Sweep Generator 

The ciixsuit for this part is shown in Figure 3 6 
The digital output (with frequency tv/ice the frequency of th 
input sine wave) of the previous stage is applied to the 
clock input of 74191 It is a binary counter with dovjn/UP 
mode control active low liOAD GNABIjB control inputs it 
gives an output pulse R at the pin no 13 which goes low for 
half of the clock period before overflow/under flow occurs 
The digital outputs of the counter is fed to the DAC so that 
we would get a digitised waveform depending on the control 
logic used to the counter For example keeping hD « high 
EN « LOW and DNAjp permanently at high or low we would get 
negative going or positive going sawtooth wave But the 
sharp fall of output corresponding to transition 1111 •* 0000 
(or the reverse for negative going sawtooth) may give rise to 
glitches or other hardware problems So we prefer to use it 
as a triangular wave generator This can be done by cont- 
rolling DNAjp mode control as shown In figure The n«3gative 
going edge of the ripple clock is used to dock the toggle 
switch (made from JKFP hy holding both J and K input at high) 
whose Q output is connected to PINS (M?/UP) of 74191 Con- 
sider the situation when DNAjp « 0 the system is counting 
up When the outputs 0^ 0^ counter attains 

111! the ripple clock goes low \dilch forces the 33N/UP 
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input to go to 1 since a - clk mx/Mn ^le m/m » 

^ 

Qji^C^C^c^ + dn/!;p The latter transition in 

um forces to go high again after a f e ? g ays Any 

way since I^/Op now latthes at high state the counter will 
start dom counting (lin * mo so on) Bafore reaching 

the state 0000 the low going transition of r will again 

c 

change the counting state v/e will thus get a triangular 
wave 

The first, one of the DACs in the system is simply 

used as digital to analog converter whilst the next one is 

basically used as an atta'^xiator to g t a wider range of bas 

steps In fact a single chip like AD390 can be used to 

attain both the functions but we have not used them because 

of their immediate non-^avail^ility to us DAC08 is a TTL 

con^ja table 8 bit DAC It works upto +10 V reference with two 

quadrant wide range multiplying capability it is highly 

linear and gives complementary current outputs ^ouT “ 

IFS for all logic states ^ere IPs * -SE (see Figure 3 6) 

REF 

These currents are sinking currents so that for taking measu- 
rement on pnp transistor Iqui. (with i^y^' pin grounded) currant 
can be directly connected to the base of the test transistor 
And to apply negative base step voltage a load resistor 
Rj^ {- 10 k) is applied between PIN 4 and ground (with pin2 
ground) and then using a voltage follower to buffer this 
voltage The range of the voltage step available would be 
from 0 to Since we are interested in testing npn 

transistors also the output is connected in the inverting 
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Fig 38(a) Astable Multivibrator (b) Waveforms across the output and 
across the capacitor, j3£:R2/(Ri^’R2)>Vo=V2+-Vd and 

R3>(Voh -Vo)/Isc 
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mode (as shovnn) so that one gets inverted positive output 
vlthin the range 0 to 'This voltage step then can be 

directly applied to the test transistor (np ^ J'e base for 
taking measureroant To get a current step from this 

voltage step we have passed it through a VCCS 

The VCCS that we have implemented here (Figure 3 7b) 
is somewhat different from the standard type (Figure 3 7a) 
Using a single opaitip for grounded load configuration The 
latter suffers from the disadvantages firstly of poor common 
mode voltage rejection; secondly the output resistance of 
the controlling voltage source Influences th^ adjustment and 
in addition the current supplied by the oontrol^voltage 
source Is dependent on the load resistance; thirdly any 
in^rfect matching of the resistances R 2 dQQtades 

the output shunt resistance of the current source appreciably 
Our circuit is more favourable in these respects 

To analyse the operation of the ciscouitry (Figure 
3 7b) we note that the second opait^ is simply acting as a 
voltage inverter So applying KVL for the circuit w3 get 


o - V, « V4 + ~ V- 
2 1 1 Ij 


The application of KCl to the output gives 


(3 7) 



Elimination of V 2 from (3 7) and (3 8) yields 



(3 8) 


(3 9 ) 
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The output current will be independent of output voltage if 
the condition 

^ ^2 (3 10 ) 

is satisfied 

In our design we have chosen ^ resis- 

tances and axe derived from a single 10 k. pot so that 
+ ^2 “ ^ * ^2 satisfied Input dc offsets are 

appropriately taken care off A small capacitor (100 pP> is 
applied to the input OPAMP (shown in figure) to avoid oscill- 
ation Note that the choice of low (and hence 1^) is 
restricted by the output currant drive requirement of the 
second opatnp Again choice of high Rj^^ (and hence low R^) may 
create excess voltage drop across it and the opanp may get 
saturated (this is also evident from (3 7)) Ihe usable 
range of R^ for a given can be worked out 

3 5 CONCLUSION 

The circuit design we have discus sea so tar is 
assembled and the system worked out well it ensures preci 
sion measurement which is so important for modelling second 
order effects For example in the characteristics 

with Ijj held constant normally shows a slight upward slojpe 
in the active region because of Early effect (vide section 
2 3) Now even quite a small variation of is sufficient 
to mask this effect Hence a constant current source is 
mandatory Mary variation of the circuit is possible 
suppose if orte wants to display single or a restricted set of 
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curves th^t can be done by proper logic conA®ctlons to ld 
and UP/DM input of the counter The only disadvantage of the 
design is that it requires a signal generate This depen- 
dence can be omitted by the siitpla additions! design as 
sketched in Figure 3 8 The output synwie tribal square wave 
as before has to be processed through ZCD in order to make 
it TTL compatible The charging and discharging voltage 
across the capacitor is first buffered and then can be used 
as input voltage to Figure 3 5 as any periodic waveform 
is sufficient and its shape is unimportant Or one can 
derive the waveform from power line itself via transformer 
or varlac ihough not meant for this system can also be 
used for monitoring other devices like mosfST tJPET etc by 
suitably connecting to appropriate outputs 
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CJiAPTER 4 

FORMUIATION OF THE MODEL 


4 1 INTRODUCTION 

Modelling t} 7 & art of c/iaractarlaatron of the 
behaviour of a physical process The devic model attempts 
to describe the tetminal electrical behaviour of the dev ce 
and the user need not know the internal phys cs There are 
many ways of approachincf this problem depending on one s 
knowledge on the internal mechanism and on the measurement 
techniques In case one has little idea about the system 
otie can propose an empirical model and the parameters of the 
model Can be obtained from experimentally measured behaviour 
of the device as viewed from its terminals This involves 
application of curve fitting techniques to obtain functional 
relationships between the terminal quantities of interest 
The result itself may help the user to get insight into the 
system But the problem in using this type of model is that 
the operating and environmental conditions are identical to 
those iictdtejr f-he merasarements have been made The 

physical model on the other hand is based on the analysis of 
the basic Internal physical mechanisms and the parameters 
chosen naturally would be related to material and structural 
properties of the device Hence the model equations would 
themselves remain Qanonical with respect to changes in opera- 
ting conditions as the exogenous variables (eg temperature) 
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are automatically taken into account in other words with 
known terminal behaviour they enable th« user to extract the 
nformation on the oper ting condition also 

A third class of techniques that we can think of 
(of course after the advent of digital computers) is r ost 
robust and begins with the basic set of equations ( like conti 
nuity equations Poisson s equat ona current equat ons etc ) 
along with suitable boundary conditions Though the solu 
tion sums to be most e^tact it values a lot of computer t me 
and memory often making it unv able and unrealistic As far 
as device people and design people are concerned this does 
not provide ny substantial benefit The model e have 
described so far is a physical model developed frwn the 
behaviour of semiconductors at the junctions and the bulk 
reg ons 

In short our approach would be like this The 
transistor operation in different well defined regions is 
considered Different approximations valid for tne respec- 
tive regions are called for to attain an lytic expressions 
These are then solved on the computer The computer result 
is monitored at various stages to check the self -consistency 
of the aPP^cximations and to check whether there is any 
hidden critical assuinptlon that m y arise the model to 
diverge from reality (validation process) These d-mand a 
particular set of experiments in a particular region of 
operation which provides bhe input data Anyway the para 
meters extracted from a typical set may be required to be 
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used to find out its associated parameters in another run 
Hence some sequences have to be maintained and one may have 
to terate amonost various stages In Apj>endi V we shall 
first discuss the least square curve fitting method (bevenburg- 
Marquard algorithm In Section 4 3 we will first summarize 
the basic equations borrowed from Chapter 2 and discuss the 
relevant ej<perimental techniques In Section 4 3 we will 
show the way to e^jploit the more sophisticated model Before 
completion we will indicate the compatibility with the 
existing SPICK simulation program 

4 2 MODEb EQUATIONS ELECTRICAL MEASUREMEt^S AND EXTRACTION 
OF MODEL PARAMETERS 

A model is of little use if acquisition of model 
parameters is difficult In this section we shall discuss 
a list of simple dc and ac measurements along with the nece- 
ssary set of equations that are sufficient for the extraction 
of the model parameters Each set of measurements tends to 
emphasise a separate physical effect and allows extraction 
generally of three to five model parameters in cases where 
a physical effect is significant in more than one set of 
experiments an iterative process is necessary For our 
case we take 2N2219 as the test transistor it is a medium 
power switching transistor and can be use^ias an eunpllfier 
also Some sets of data have been taken straight away from 
the data book [5] for the verification of the software 
developed 
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4 2 1 aa a Function of and as a Funct on of 

As has been pointed out at the beginning of SectJLon 
2 3 1 the usual estpression of junction capacitance given by 
(2 19) suffers from singularity This has been bypassed in 
the proposition of ( 2 19a) So if one takes its derivative 
the following result is evident 


oe 


( + b) 




( 1 + 


n. 


^ “ *^e + b 


(4 1) 


where x The elements of the vector P in 

(2 19) is then given by 

Pi ° ''bib P4 = P2 “ V® P3 “ 

@ 

(4 la) 

The initial choice of the parameters are important A rough 

estimate for them can be obtained as follows One may use 
N N 

Vfijg « Vy In or it can be assigned the default value 

"i 

0 7 V in case the doping levels are not know (cel ed 

the grading factor) is estimated from the slope of 

In Cq versus ln( 1 - ^oe measured sero 

bias capacitance A forward biased capacitance can be 

estimated from the extrapolation of the slope 

of the reciprocal cut off frequency such that 

d( l/f^) 

°ef “ 2H dTv^ 

a ^2 related to by [39 ] 


U 2) 
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b 





oe 



(4 2a) 


where r is constant close to unity Itsf exact v lue depends 
on doping profile These estimated valu^can be used aS 
initial guesses for a non-linear curve fitting procedure 
which fits the measured values of to the calculated values 
by adjusting n^ and according to (4 1) Note 

that the parameters finally obtained are purely model para- 
meters though they can be attributed to jfhysical definitions 
as discussed above 

By a similar procedure one can cS^-culate n^ 

C and and in turn the elements of Pq before 

that one thing is to be noted The terminal capa- 
citance is the total capacitance ^i^ whereas 

represents the capacitance of the active base- 

collector junction Anyway all the re^*®tions in (4 la) can 
be used excepting the third one which six^wld modified to be 


^oo^BIC 1_ , (4 lb) 

P3 “ 1 - 1 + 

where r^ is defined in ( 2 l23a) 

A 

4 2 2 hpg Ver us (At Low Injection) 

In Section 2 2 we have deduced tr» basic equations 
relating junction voltages with currents m the normal 
active mode equations (2 10 ) and (2 13) j^’educes to 


^ i S’ 
^cc 


Ig 


(4 3) 
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and 

- V^p +isE* ' ' 

where currents and voltages are in normalised form (eg 

'^be 

ic *= '^be SUmlnating e from 

them we get the following useful expression 







i/n_ 

) Vi. 


where 


-i/n 


SE 


^SB^S 


e 


(4 5) 

C4 5a) 


This celablonshlp can be utlllseci to extract the 
parameters I,^ Igj. and from hpg versus curve pro- 

vided we have the knowledge on before hand For this we 
take the equation ( 2 9) divide throughout by and assuming 
Q B T^Ii. (by definition of I,^) we end up with the following 

expression for gj^ 

% = 1 + «be + '*bc + ®^c 

where we have neglected the last term The expression for 
^be '*bc “vail^le in ( 2 20) and ( 2 21) The rele- 
vant parameters (Pg P^) has already been extracted from 
reverse capacitance voltage curve (see Section 421) The 
only unknown quantity left is B By carefully choosing 
(I V ) points such that remains less than the critical 

current (funoUon of see Section 2 3) at which base 

push out effect starts occurring enables us to put B - 1 
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After getting the above four parameters we now 
concentrate on the tail porfion of the hpg curve where 
the base push out effect ia predominant By Introducing 
appropriate expression for the f actoc^B in ( 4 6) the charac- 
teristic parameters are then amenable for numerical evalua 
tion But the push out effect also comes into Play in the 
experiments (4 2 3) and (4 2 5) This means that data from 
these measurements have to be fitted alternatively until a 
consistant set of model parameter is obtained we shall 
describe the procedure in Section 4 3 1 

4 2 3 Measurement of Emitter and Collector Series 
Resistances 

Shackley et al [42] have Indicated that for an ideal 

transistor in the common base configuration the f lo ting 

(I =* O) collector to base voltage will be almost identical 

o 

to the forward biased emitter to-base voltage But In an 

actual translator the difference between the two voltages Is 

attributed to the volume recombination of carrier pairs 

which diffuse from emitter to the collector Accord ngljf 

In an Ideal transistor In a common emitter orientation 

V t = O whereas in an actual transistor with emitter 

11 » O 

series resistances Rg « o ** ” B 

larjy If the emitter and collector are interchanged to 
operate in eotmnon-collector orientation then with a collec. 

tor series resistance Rq + ■ o * ” 
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Thus measurement of floating collector and emitter 
voltages as a function of forward biased base current will 
provide evaluation of Rg and These measurements h ve 

been done by our curve tracer as follows The base current 
step is as usual applied to provide the forwcord biasing 
whilst emitter is grounded and the collector is kept floating 
Then the collector and emitter (ground) points are connected 
to the measurement unit ( see Section 3 3 2) to measure the 
voltage (V^g) This ensures a high impedance («^1 GO) 
essential for such measurement For measurement of ^CC 
above connection has to be made in reverse Then from the 
display the slope of the lines would give Rg and 
4 2 4 as a Function of V^g at Constant V^g 

If the convention for terminal current that all 
current flowing in is positive Is adopted then the terminal 
BE voltage is related with true junction voltage by 

VbE “ ^BE V'Hj + 

Where Ig « -dc + emitter current The symbols 

have their usual meanings Cortbining (4 7) with (4 3) we get 


Vbb ' '"t + V-’b + 


(4 8) 


If is not measured along with Ig then one can use 
instead of Ij, in (4 8) where h^j. le to be taken either from 
previous experiment or through the use of (4 5) since para- 


meters lik& 


fl 1 are already known Bguation (4 8) 
Ff ^SE 
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can be subjected to curve fitting technique to squeeze the 
parameters ig and is very small normally of 

the order of 1 o (the contact lead resistance) Mote that 
we have not accepted a standard value of C 25 9 mV at 
300 K) Since the operating temperature is unknown and it 
appears in the exponent it has got significant effect on the 
magnitude of the current and thus there is enough justifica- 
tion to take it as a model parameter At low injection 
G 1 and it is not possible to get and r^^ separately 
(because of lack of sensitivity) rather their sura » 

-j. should be taken as a single model parameter At high 
currents where the base resistivity modulation is s gnif - 
cant one can use the same equation (4 8) to extract rj^ 
alone accepting the already known above parameters 

4 2 5 Output Chataoteeistics (V^g versus l^,) in Saturation 
Region 

Here the BB as well as CB junction are forward 
biased Hence both of them have significant r 1 e e on 
the currents in the system For the time being we are neg- 
lecting the oontrlbution from overlap diode and conductivity 
modulation in the collector region we consider the following 
set of equations(4 7) ( 2 13) ( 2 64) ( 2 10) and the two 

given below 


- v„ <is. * ■ ° “ 

oan be solved foe the other four 
T V these equations can ne ew* «« 



96 


variables All the constants barring \ 

already known Using their known values and accepting some 
initial value for the remaining four one can solve the above 
set of equation in a manner already prescribed The calcu 
lated value of is then compared with the measured value 
until beat fit is obtained The parameters are updated in 
each iteration until the convergence is reached Alterna- 
tively one can extract the same four parameters by going to 
inverse active region and then adopting a procedure similar 
to that written in Section 4 3 2 
4 2 6 f^ as a Function of at Constant 

This experiment can provide a lot of information 
for high frequency and transient operation of a transistor 
and also it helps to determine the modal parameter of base 
push effect But there exists a lot of confusion [38 39] 

about the exact definition of f,j. and the way it has to be 
measured Hence a cleat understanding of f^ is necessary 
liSt us define the transit freguenoy ^ 

f » ipif 11) 


where p = short circuit common emitter current gain 

From this it is to be seen that if IP I varies inver 
sely with frequency (or it falls with 6 dB per octave) then 
f is independent of frequency and it equals to f^^ the 
ulity gain frequency if this approximation is vaUd then 
it is easy to measure f,, measuring p at any high frequency 
and then taking help of (2 U) But we shall see below that 



the permissible frequencies for this measurement ate restri 
cbed to 3 cortsdn 

The conventional expression for the low frequency 
dependence of p 

where i; is the emitter to collector delay time writing 
1/fp for 2 wt the modulus form of (4 12) appears as 

•pi - Pq /l + (f/£p)^ (4 12a) 

Here fp physically stands for 3 dB cut-off frequency Defi 
ning the ideal transit frequency by 

fTo = Po*p 

We get the following relation to from (4 11) 

f^ = 14 14) 

/l + (£^q/£ P^)^ 

We can at once conclude about the lower limit of the usable 

2 

frequency range Assuming (Ipl/P^) <<1 we see from 

(4 12) and (4 13) = (fg/£)^« 1 and equation 

(4 14) reduces to f^ « f^ Thus the measuring frequency 

chosen should be at least some five to six times larger than 

fp The upper limit to it (over which f^ wxll again become 

frequency dependent) is caused by parasitic transistor and 

pactcaoe parameters i e they contribute to higher order 

lower 

poles and zeros which are far^than the cut-off frequency f^^ 



itself This limit is esloulatsd in a laborious way in the 
paper [403 and also the optimum frequency is suggested 

Let us now turn out attention how to extract infer 
raation from the f^ versus curve one can define the 
emitter to collector transit tiit» 


X 



<4 15) 


where 0^ represents the total excess majority charge in the 
emitter base and collector Junctions and quasl^neutral regions* 
Using quasi -neutral condition one can also equate excess 
majority carrier charge with the excess minority carr er 
charge In the active base region the model charge Is 
And its derivative with respect to gives the following 
terms 


<30b 

+ ^^c + «E + ‘‘c] + ® ■'£ ‘“ 16 ) 

where we have taken help from (4 3) (4 9) and (4 10) 

^b " ^1 **' ^2 other symbols are already defined 

{Actually to be replaced by x where lA^j^ « 
where is the recombination life time in the base) i^ote 
that this is valid at both low and high injection If we 
take into account the capacitive effect of the inactive base 
region it provides an additional delay 


fVV 

«c' 





t4 17) 
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where is definejin ( 2 123) and is gx en by 

Vc 

then (4 17) can be evaluated to be 



Vt 


''ife +«c ^•“e +«b^ 


(4 18} 


(4 19) 


If we sum up (4 16) and (4 19) and also include n that the 
parasitic capacitance Cp^ between the emitter base terminals 
(it consists of the capacitance of the bond pads of the pack- 
age and possibly of the measuring socket) which acts n 

parallel to C we end up with 
6 


T « 


2n £T. 


« B T 


fN 




C4 20) 

Here C* » (4 20e> 

+ °pe^ + ^W% 3/ P > < 4 20b) 

The last term in the expression (4 20) incorporates the 
delay time in CB junction depletion layer w la the deple 
tlon width given by ( 2 24) Is the saturation velocity 
Usually this term is very small and can be neglected 

Let us try to explain Figure 4 1 in the light of 
( 4 20) Ths shape of the curves reports three physical 
effects At low collector current is dottlinated by the 
second term 4.n (4 20) which is ^proximateiy C At 
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low collector current the difference in f^ between the 
cases with = 5 v and with « 10 v is resulted mainly 
from the first term in (4 20b) The maxinum f^ is limited by 
which is roughly Wj^ 2T1 At high collector currents 
base push out effect causes rapid decrease in Hence at 

low but different dataon f^ helps in determining 
since + C^) is already known from (4 1) at inter 

mediate values of an estimate of is possible Fran 
the f„, (large) data the model parameters of base push out 
effect can be found out 

AS already mentioned at the beginning of this sub- 
section a necessary condition for determining the para- 
meters from the measured curve ® Ipl f (equation 

(4 11)) in the manner described above is with 

from (4 20) That means f^, must not depend on the measure- 
njSit frequency within a certain current range First of all 
we note that for a diffusion transistor the transport factor 
«= Sech (Wjj/Iijj) the frequency dependence of o(w) is 
therefore given by a(w) Sech i” 

base diffusion length is replaced by its complex value ■= 
Lj/(1 + jw In view of the relation p{w) = a(w)/ 

(1 - tf(w)) we see that p(w) is no longer given by ( 4 12) 
rather it is approximately valid if one expands the Sech 
term and retains the first two in the series when base 
grading is present an empirical expression le proposed. 

«o jmw^tt <4 25) 

a(w) = i + jw/wg 
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where w <=» 2itf 


Is the operating angular frequency 


with 


and 


ra « O 22 + 0 It) 


wa 


0 - + (n/2)^'^^] 


(4 25a) 
(4 25b) 


At low w the exponential term in (4 25) reduces to unity and 
again (4 12) can be realised At high frequencies the origin 
for the strong dependence also comes from parasitic lead 
inductances and stray capacitances Figure 4 1 clearly 
displays this feature it is also obvious in the low freq- 
uency range is nearly independent of frequency and 1/2 iif^ 
Increases almost linearly with l/l^ in accordance with (4 20) 
At high frequencies however increases strongly with f 
The result* for ~ 50 and oo are Included in the 

same figure Since the frequency chosen is low enough to 
avoid high frequency effect the curve with 6 « cd shows 

^To "View of (4 14)) and hence an ideal one We 
therefore note that the best choice of frequency range is 0 3 
to 0 5 GH* for this high frequency transistor We also note 
that since high frequency effects start becoming significant 
at much lower frequencies than itself (i e the presence 
of secondary poles at much lower value than actual 

unity gain cut-off frequency significantly differs frcjm (4 20) 
and also the phase lag is considerably different 

So far w® have not included the storage charge 
effect in the emitter and collector regions In the emitter 
region it creates a delay ^ where is the 



excess minority carrier and 1^(0) is the hole current in the 

emitter and represents the ratio of the effective oummel 

number in the B ami b r-^ 4 ^ 

B regions respectively The expre- 
ssion is evaluated and found out to be 

® Vn “^othUgA ) ■ (4 26) 


The collector storage effect is important in satura- 
tion and inverse active mode operation We have seen in 

( 2 93) that Charge dens ty falls Unearly with distance so 
that ^ _£1B „_j ... 


“ H “^^"9 this approximation in the third term 
of ( 2 86) multiplying throughout by qA^, and then integrating 
over the limit x s o to x w^ib which conductivity 


modulation is important) we get 

PcE I qAc(p(0) + Mep)W(,iB 


(4 27) 


taking the approstimation p( 0 ) » and utlUsing the 
relationship <2 87) wa get an appro3<lraate form of the coll- 
ector storage delay 


T 

a 





(4 28) 


where is to be chosen from (2 95) These two delays 

( ^nd % ) have to be Included in (4 20) 
e c 

One fact has ltd seed our consideration so far we 
have neglected the last term since Vg is large But at high 
current densities the electric field in the depletion region 
reduces because of Kirk effect and appropriate expression 
for carrier velocity ought to be used from (2 6^) to ( 2 68) 



104 


instead of using Vg only Hence this term gives rise to some 
contribution to the total delay at high collector currents 
Before passing we should mention that the softwauce 
that has been developed is not based on the closed form e^tpre- 
ssion (4 20) Rather the set of sl 5 « equations mentioned at 
the beginning of Section 4 3 5 along with (4 18) where 0^ = 

Ob + 0^ picked up They include eight equations in 10 
variables Q^) 

Given and they are solved for by N R techniques for 
the remaining eight variables Furthermore the first deri- 
vative in (4 15) can be conputed By adjusting the value of 
the model parameters the calculated value of is fitted to 
measured values The parameters extracted will be discussed 
in the next section 

4 3 MODBLliING ON THE BASIS OF EXTENDED GP MODEL 

In the previous section we have discussed the method 
of measurement and extraction of pure GP parameters In this 
section we shall take the help of the same set or measure 
ments with one or two new additions to compute the relevant 
parameters We have prescribed in Section 4 3 I to make 
Junction capacitance-voltage measurements and thereby to 
extract the parameters (P^ P^,) which ate required to model 

Early effect We shall sat aside these parameters and 
Instead include Vp^p Vpgj, and Vp^^ defined in equation ( 2 26) 
( 2 32) and ( 2 37a) respectively as the model parameters for 
the same Consequently the set of measurement one have to 
maice radically gets changed 
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The measurement of Is done indirectly a slow 

ramp is applied to the reverse CB junction while its capaci- 
tance is monitored At voltage v and more than V the 

i “ii' rep 

reading should show a constant value because the depletion 
region in the collector side touches heavily doped (n’*') 
substrate and to the base side it is going along upgrading 
Hence the onset voltage for constant C should be identified 
with This experiment is performed using PAR410 ramp 

generator with HP7015B XY Plotter for recording 

Measurement of is directly possible if this 

voltage is less than CB junction breakdown voltage 
In that what have to do is to short the BB Junction and 
Increase (reverse) until the base depletion edge touches 
the EB junction The CE current will show a steep rise at 
this voltage (Vp^p) which can be monitored In case it la 
not possible one then can estimate it from the slope of 
I^<0)/I^ versus + ^biC where l^(o) refers to coll- 
ector current at « 0 The parameter s not so much 

imgportant and One can follow the same procedure as is used 
for the determination of Note that the determination 

of Vpjp V^^p allow the determination of N^p and w^p 

They can be found from CB junction measurement also but 
the parasitic effects impair accuracy 

one point we want to elaborate here is that the 
ejtpression for that vouW be obtained from < 2 33) 

or ( 2 34) ae it suits should be Inserted in ( 2 42) to get 
effedtive normalise charges associated with CB Junction 
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wllih slight change we rewirlte this e^tpresslon here 


^bc ef£ 


N a 

^ v) (1 + a)( 1 

1 ■> 


I + (- 
1 




+ q 


a 

be 


(1 + c) 


(4 29) 


where what we need to know first hand 

±s the values of n^; other i|uantitles are already known 
They can be estimated from device structure We note that 
too much accuracy in their values is not needed here as the 
modification (4 29) itself has little sign flcance 

The ne3<t set of parameters Py Ig^ and remains 
essentially unchanged only the related equations are somewhat 
different Note that in (2 104) we have introduced the quan- 
tity and where ^2 " ^SB write 

i l^i (o)K where K is to be found out from ( 2 lo9) or 
cc r c r r 

( 2 119) then the following equation can be derived 


1 

“fb 



F n»l 


^n 


(1 (o)qv,(o)) ® 

i^CoJiT 


n«l 


en 


(4 30) 


which is similar to (4 5) The way of deducing the values 

of K k s etc have already been described in section 2 6 

r n 

we also want to specify that as a whole has to be 

calculated either from ( 2 63) or from ( 2 65) In them the 
part q^^ would be obtained from (4 29) and (2 37) while the 
part by first determining (2 62) via (2 61) 

and ( 2 60) and then Putting in ( 2 63) or in (2 65) in t-he 
latter case The third way in tackUng the problem which 
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we feel more convenient Is to adopt lumped parameter model 
(jpagure 2 13b) Here we take the simpler form (4 5) as our 
objective function The base series resistances taken would 
be and R^ the latter one is to be determined from (2 113) 
or { 2 114) or from (2 120} as the situation entails In 
the part would be obtained from the prescript on sugges 
ted while dealing with the equation ( 2 122) or from ( 2 63) 

But in the latter case an average value of has to be 
substituted 

The third set of parameters re Ig Rj^ and 
For their evaluation the technique dlsclos d in Section 4 2 4 
xe recommended 

The fourth set of parameters nameljf Ig^ 

T ( « tVt ) would be identical to those in Section 4 2 5 
t r r 

if we take the definition of as in ( 2 132) where the 
aontribution of overlap diode Dj^ automatically enters But 
if we accept the more sophiatloated model (Figure 2 13b) 
then the presence of R,, separating the intrinsic part of the 
transistor from the overlap diode has to be considered The 
definition of would then come on an equal footing with Pp 
i e the term Containing 4 has to be deleted from { 2 132) 
The set of equations those which will be relevant are ( 4 7) 

( 2 13) ( 2 10 ) ( 2 126) (4 9) (4 10 ) (2 130) and (4 18) 

But because of the presence of the current Ip^ the following 
modifications has to be made to the right hand side of these 
equations before we can use them, ta) subatract a term 

from (4 7 ) (b) add and subtract the current Ip^ to 
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(2 13) and (2 10) respectively (c) add to 

(4 10) (d) replace by in (4 10) (e) subtract 


to the last of these equation the rest remaining unchanged 
Note that we have chosen equation (2 130) instead of (2 64) 
This set consists of eight equations with 10 unknowns The 
concerned parameters can be evaluated following the procedure 
adopted in Section 425 

4 3 1 E:<traction of Parameters on Base Push Out Effect 
This requires separate attention as this ffect 


appears in more than one set of measurements and the informa- 
tion laid in should be exploited for the simultaneous extrac 
tion of the parameters We have mentioned in Ghaptes 2 Sec- 
tion 2 4 that the factor B in (2 63) should be replaced by 
B = appropriate expression 

for w-,n available from the set (2 75) - (2 77) as the 
situation demands By introducing the terms - ”^p/Wb 

is ” ‘^'^o^ep^S '^ep ep S 

results are derived from the above set of equations* 


(a) ^eJg 


B 




tr J 


ep 


, Js ’'BIC -BS i 


J- 2 V— 2 


{4 31) 


B 


(4 32) 
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(c) 


•^C ^ 




2(V 


BIC 




^2 
s "s 


(jc 


(4 33) 


Thus the parameters chosen are r t j_ and V_„ Their 

w 6 5 ep 

initial values can be chosen from a knowledge of device 
structure and material properties the relationships being 
given in the first paragraph of this section Since the 
factor B appears in both (4 6) and (4 20) the parameters 
evolved from one set of measurements re used as i itial 
choice to the other set and this is done alternatively until 


convergence from either set is ensured 
4 4 COMPARISON WITH SPICE MODEL [ 35 ] 

We have included this section with the view to 


make it (the proposed model) more convenient to the users 
for people ate generally more familiar w th SPICB model 
in the modified form of SPICE aa much as forty parameters 
are ehosen They are Usted in Appendl:< III There a 
separate column has been added to show the compatibility 
with the present proposition Before going for any compa- 
rison we Should mention that in the software package the 

Til etc are taken as TP TT 
name of the parameters \ 

IS, the reason is obvious Those which ar identical are 
indicated by + a gh Those which are not same bat can be 
a derivative of our model parameters the relationships are 
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indicated Again SPICE program takes care of some more 

phenomena which we have not taken As for example it has 

included the temperature dependence of I- via the terms XTI 

5 

and EO; the presence of junction capacit nee (between the 
collector epi layer and the substrate) via CJS VJS MJS etc 
They can also be incorporated straightway in the present 
model It is found convenient to use expression for some 
variables instead of modelling them by more than one para- 
meters For example the variation in forward transit time 

has been taken care by introducing three additional 
r 

parameters XTF VTF ITP in SPICE However we are making 

use of an expression { 2 62) for th reby avoiding the 

necessity of additional parameters Hence fewer number of 
parameters are required in our version Apart from this 
we have included some more effects like emitter edge 
crowding mobility variation in the base region base push 
out conductivity modulation in low doped collector region 
etc which have been sirrply ignored in the modified version 
of SPICE It is iso to be noted that ours approach in 
tackling particular effects considerably differs (for 
example SPICE has accepted Hauser s expression [23] for 
modeXUng the modulation of base resistivity but we shall 
incline to use (2 114) or (2 120) for the same) 
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SOFTWARE DEVELOPMENT 


5 1 Introduction 

In this chapter we will first briefly described 
the parent equations that are taken s objective functions 
of our optimisation problem These equ tions are in fact 
already discussed in Chapter 2 and Chapter 4 but in a 

different forms The reformulation of them is done as well 
as a different set of parameters is chosen with a view to 
attain greater numerical sensitivity The complete flew 
chart and the program itself will be treated separately 
in Appendi3< IV In section 5 2 e 3 <perlmental graE*i and 
data will be discussed vhereas xtracted output study 
will be the subject matter of the next section In 
conclusion we shall indicate about the convergent 
criteria the reasons for probable failures of the program 
in pathological cases and the way for their bypassing 
and remedies 


5 2 Qbieotiva Functio ns for the simulaUgn 

The listing of the parameters of the. simplified 
version o£ the model aret ifbe eight parameters (Pg Pj,) 
related to junction depletion layer charges; saturation 
current space charge currents Igg IgC forward and 



reverse ideal amplification factors p j, and p the 
forward transit time u „ transit time ratio t ^ 

(» and forward knee current required to model 

e5<cess minority carrier charge in the base region; the 
four parametets r^ and v^p related with base 

push out effect Apart from them some more parameter s 
resulting from the parasitica are also considered to get 
the terminal behaviours of the transistor 

The inport nt ones are and resi tances 

of the active and inactive past of th base; R_ 

resistance due to collector epilayer and that of the n*^ 

substrate; the parasitic capacit nee C_ between base 

pe 

and emitter contacts (as discussed in section 4 2 6) 

r* the ratio of the active to inactive part of the collector 

A 

junction; J the saturation currant of the dio<i- D, (see 
Pig 2 13(b)) 

The four parameters are related to 
and 3^2 via equation (4 2) wher as the letters govern the 
Junction capacitance by (4 1) The way of extraction 

is indicated in that sectlpn (i e sec 4 2 1) But from 
the optimisation solution point of view the equation is 
highly eccentric and crucked and the idea of extraction 
of four parameters may be sometimes misleading Any way 
we take the logarithm of the expression and use the new 
set of parameters defined below so that the equation beccmies 
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of less crucked nature 

= P(3) - P(2) log . log (1 + 5 -^, /iil.) 

B£*N 

C5 1) 

where - 1 + P(l) In C^; the subscripted 

variable P Is simulation parameter and not to be confused 
with Pg listed above tts elements are related as 
P( 1) = lAgjg; P( 2) n^/2 P( 3) In and P( 4) « 3^2 

In general fitting a set of points v^g) onto the 

theoretical expression would enable one extract the 
parameters P we note from (5 1) that the expression is 
rather insensitive to the change of the value of P( 4) 
this may result in ill-conditioning of the coefficient 
matrix and extraction of parameter would generally beconve 
d fficult To get rid of this problem we recall that 
the Parameter P(4) comes in our way as a mean to avoid 
singularity in (2 19) Physically this models the presence 
of finite capacitance in forward bias when Vgg 
Hence if sufficient data in forward balsing conditions are 
taken the insensitivity can be minimised 

A similar description follows for the extraction 
of the elements of P_ But here the value of P(4) is 

C 

very small and it is advisable to replace P( 4) by the 
constant value a o which can be estimated fr^n an 

C £ 

expression like (4 2) 
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The third set of parameters (V_ R. . ) re 

T o OD 

identical to that listed in sec 4 2 4 the parent equat on 

being 4 8 For this we have to have knowledge on and 

/ 

hyg The appropriate values of the latter can be picked 
up from (hj.g l^) measurenvent For we t ke Its expression 
(4 6) ; in that q|^ h s to be found from ( 2 20) by 

using already known parameter from the last two runs In 
the fourth term we put B 1 nd for the value of we 
have to rely on experimental curve It is advis ble for 
extracting this set to t ke data points from low V^g 
portion of the curve Fcr thi choice inv Udat s the 
approximation B»1 makes the magnitude c£ the fourth 
term in (4 6) insignificant so that an estimation of 
the value of Ij^ would not lead to any ser ous error and 

also the total base resist nce ^ 

constant and a maxiraum since base reslstlv ty modulation 
is absent On the other hand too much low value of 
current should not be used as then the second t.rm in 
(4 8) become Inslgniaoant nd this may lead to difficul- 
ties in numerical analysis Hence a compromise has to be 
made in choosing the data points 


Next we turn our attention to sec 4 2 2 
reasons already mentioned, we reformulated C4 5) 


For 

as 


(5 2 ) 
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and 

% *®bc + 

where the simulation parameters P s are P( 1) - 1/pp 

P(2) « Igj. P(3) = 1 /t)^; P(4) « 1/lj^ the parameter 

I is already e 3 <tracted in the previous simulation The 
s 

parameters can be e3<tract d as usual Note that it is 
not hpg but the reciprocal valu has to be used as in 
put to the program 

The measurement on Rg and is directly possible 
( see sec 4 2 3) and needs no simulation 

The parameters p Ig^ (important 

while the CB junction is forward biass =d) are discu sed 
in detail in section 4 2 5 for their simulation procedure 
The parameters chosen for our programme are PC 1) - 
P( 2) « I ^(3) * Vr^ reason is obvious 

Their initial values can be estimated from their expressions 
derived in Chapter 2 

The last set of parameters is on base push out 
effect evolved from simultaneous fitting then to (5 2) 
and (4 20) at high values of oollectot currents TOUs 
matter has been treated separately n section 2 4 2 and 
section 4 3 1 which is sufficient for the understanding 
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of the program By estperlence we have seen accepting value 
of from material properties rather than to 

use it as model parameter gives better numerical sensitivity 

5 3 Results and Discussions 

our result and discussion would be based on the 
sliAplified version of the model yet enough accurate. The 
software developed is first tested by talcing points from 
characteristic graphs supplied in the data book and then 
from measured values obtained during the course of the 
project work The end results are seen matching; some 
differences is also observed the reason for which is also 
predictable 

In the page 117 of this book is shewn a reprint of 
the relevant graphs of 2N2219A siUcon NSM transistor from 
tfatronal Semiconduetor Data Book [ 5 ] Out of them we 

have only used Fig 5 1(a) (b) (9) ^ 

fact faced a lot of problem while using them Since the 
graphs are conjusted the reading taken from them is not 
enough accurate and this leads to numecloal problems It 
is seen while optimising the objective function the value 
obtained from simulation is matching with the input data to 
a degree higher than the degree of accuracy of the reading 
yet the program is not converging Secondly the set of graphs 
available Is Insufficient to record all the parameters As 
far as ex*4>le graPh 5 1(h) should he extended enough to give 
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an idea on Knee current and knee voltage which are so important 
for setting initial choices in the program The graph of 
output oharaoterl sties at saturation and at active region is 
also missing from the list we have still used that 

because of its ready availability 

In the following table we will show data obtained from 
measurements In Table 5 1 is shown data point extracted 
from SB and CB Junction capacitance measurement This is 
done by applying a ramp (slow enough to produce any charging 
and discharging current through the capacitance) super 
imposed by a small alternating voltage (»« 15 mV at 1 MKz) 

For this model PAR 410 ramp generator is used To monitor 

the reading an XY plotter is connected which straightway 

plot the required gral* A digital voltmeter reads the voltages 

Prom the set of values of the extracted parameters 
we note a significant variation in their values for various 
set of measurements Ohough convergence is ch±..v d in all 
oases the result shows the importanoa in choosing the region 
of data points In Table 5 l( ) and Table 5 1 ( ) the 

data points are taken from forward bias region also and 
the result is more handsome and not contaminated Before 
Passing we would like to cwimenb that the wodslUng of 
■any effect through Punch through voltages etc is more 
realistic and as it consumes leas computer time it is 
aogeptable for circuit simulation 



CV Measurement for EB and CB Junctions 
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The nea<t table abovis points taten from both the 
data sheet and that through measurement Along with the 
values we have also indicated the room temperature over 
which the measurements are taken It is important to 
note the variation of Vgg with temperature at same collector 
current The corresponding table for extracted parameters 
are also Indicated 


Table 5 3 

Capacitance Parameters Values 


Parameter 

name 

Table 5 1(a) 

1 

’ Table 5 1(b) 

1- 

' Table 5 l(o) 

■ 

“■r 

Table 

5 1(d) 

. 1 

V 

BIE 

8398B4O0 

2005B400 

6i42&»00 

3652E4O0 


3757B4O0 

2289E-}O0 

2850S400 

3051B-400 


2093E4O2 

9639B+01 

1963B+02 

7995B+01 

®e2 

1130B400 

1681B-.05 

3436B+00 

64&8B^ 



CB junction, vertical scale,! div-1pF,the dashed one for 
EB junction, vertical scale 1div=2 5pF, while the horizontal 

scale IS in volts 




CHAPTER 6 


CONCLUSION AND SUGGESTIONS FOR FURTHER WORK 


In this work a unified approach for modelling 
bipolar junction transistor of either type from large 
signal dc analysis has been presented The basic 
premises of all the model expressions is from Gummel s 
charge control theory All most all known major non- 
linear effects have been incorporated nd appropriate 
reformulation is done when one or more secondary effects 
are present in particular region of operation of a transis- 
tor AS for example and high injection both conductivity 
modulation and emitt r edge crowding are present one 
Inducing the other and they have been accommodated 
collectively We have indicated the way to choose 
tacitly the data points for the extraction of a 
particular set of parameters and the initi 1 choice of 
parameters This is e3«tremely important from simulation 
point of view and as also the expressions are quite 
nonlinear in nature The s tuation in deep saturation 
and inverse active mode operation is generally more 
complex and the assumptions to he taken as a basis for 
modelling this regxon seems to be more stringent than 
that in forward active and saturation mode We have taken 
enough care for extrapolation of the previous set of expressions. 


The validation the proposed nuDdel (extended GP) 
has not been done in th^ present work vje have t sted only 
the simplified version The required experimental set up 
and the way to take various measurements has been proposed 
By taking the appropriate measurements the software of the 
program could be developed Some of the eff'Kst which we 
have not included can he easily added as for x mple the 
presence of collector epi region to substrate c pacitance 
temperature dependence of Ig etc Ther are many more facets 
of the device that one may be interested in but it is hoped 
that the above suggestions if implemented is su ficicnt 
for any practical purpose The software pack ge so developed 
may be used as a module for sinul ting circuit behaviours 
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APPENDIX-I 


BAIANCING TECHNIQUES FOR OFFSETS OF OPAMPS 

An ideal operational amplifier is perfectly balanced 
hat is whan s* 0 A real operational amplifier 

5 <hibits an unbalance caused by a mismatch of the input 
ransistors also because of temperature grading effect [43] 
his mismabh causes unequal bias currents (ig and ig 
espectively) flowing through the input terminals and also 
n input offset voltage which must be applied between the 
nput terminals to balance the amplifier 

A modified form of the equivalent circuit is shown 
n figure ( la)* Here we have taken instead of ig and ig 
he input bias current ig " input o££set 

urrent ij^^ « ig - ig when = 0 

Consider the case when the OPAMP is active in 
Inverting mode With positive input tecminai grounded 
here exists a virtual short circuit between the input 
erminals hence whole o£ the current ig flows through 
creating a voltage drop ig Rf exists 

>nlythe do offset voltage «ith posiUve terminal 

irounded appears across B resulting in a current 

he same current flows through R, Hence the output voltage 

/ = (R + H,) V, APPlyina super poslUon we see that 

O f lO- 

he neb output voltage appears to to 




g 1(a) dc equivalent circuit of an OPAMP 
(b) Universal offset balancing 
technique m inverting mode [3] 



A II 


+ ^io ^ ^ 

In case offset cureents ig anc3 i“ are either known 
or measured the contribution due to offset current can be 
ftxlly annul ie<l by placing a resistor R in the positive 
terminal so that Ig R - l" R|| (R = rI (Rj) But if they 
are not known then choose R « R so that at least the 
effect due to bias current ig is nullified leaving behind 
a small output error offset voltage effect 

Is off sett ed by using the potentiometer as shown In figure 
1(b) It allows a small voltage effectively In series with the 
noninverting amplifier in the range + V R + R ) 

Choice of R is aa already discussed in connection with 
reducing offset current effect Since Rj^ R the offset 
voltage circuitry is atiyply isolated and thus the two effects 
are controlled incJepenaently The balaneing technique can 
similarly be considered for non*^nverting mode operation 
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APPENDIX ■II 


CB JUNCTION WIDTH AS A FUNCTION OF VOLTAOB 


We Start with the one dlnvensional Poisson equation 



dE P ( x) 
dx “ e 


I p(x)-n(x)-tN^j(x)-N^(jt)] (n 1) 


and stpply it over the depletion region to the e ther side 
of the CB junction To the base side N-Cx) - N.Cx) - 

^ X * 

^ep ^ ^ collector side it is only If 

we assume the mobile carriers of the tr nsported current 
moves with the saturation velocity v^ in the depletion region 
due to the presence of high electric field then tl*eir 
number (= j/q = NgplJ/Jg) “ N^pJ^n 
defined in (4 3 1) ought to be added to the base side and 
to be substr acted from the collector side Neglecting the 
majority carrier contribution in the depletion region we 
get the equations for the respective regions as 


dS 

dx 


dE ^ 
dx 


qN 




(e 


-1 + Jcn> 51 




( 1 - 3cn) 




j 


cn 


e ’'on' 

(II 2) 

( 1 I 1 . 3 ) 


where = l-Jcn 

is B = Bjg at X = 0 


The boundary values fot the elec±rie field 
B o at X » B « Bgg at X - where 


the depletion vrd.dt.h to the base and collector 
side respectively Using these boundary values ( ll 2) and 
( II 3) can be solvedto yield 

(a) For -X < 

P 

E(x) « - 

(b) For 0 jC X < x^ 

qN^_ 

E(«) - (« - *„) 

where 

B « J v^/u 
cz ’^cn s' 


(II 5) 


(II 6) 


X < 0 

qN 


SE_ [ e * ^ - e 


•p 


- ion ^ 

(II 4) 


The maxinuam electric field is given by 

c. 


^ = “ 




e 


.1 


.) + C V 




e n oz 


(II 7) 


This gives the first relationship between and x^ We 
note that the linking current reduces the maximum electric 
fxeld Equality ( II 7) is an approximate one because of 
presence of \ <3toQS not exactly occur at x « 0 Any 

way for low current (when no base push out effect) (ll 7) 
is sufficiently accurate 



A-V 


The total voltage drop •• E( x) d>;; 

If we break up the integration limits from - Xp to 0 0 to 

then to yi and use expressions (II 4) (II 5) and (11 7) 
n ep 

in the respective regions we will land down with 


'^BIC ” ^BC ” "* ’*■ 2 " 2^cn "" 


®cz '^ep 


P 

(II 8) 


The charge conservation across the junction gives us 


- 1 


/ 

-X 


n 


<1 - 


n (x> 


e 




-)dx 


(II 9) 


P 


ep 


where ng(3<) would be given by (for where r^ is 

defined in Section (4 3 1) 


ne(3t) 


(ip E(5()+Vy dl/dx (ln(x)) 


(II 10 ) 


where n ( x) is to be solved by NR technique with boundary 
e 

values as at X - - Xp 

charge injected from the BB Junction and E(Xp) 

jlt the other end ng(x) «■ N^p and B (x^^) « - Now the 

built in electric field Vgjc ia given by 

Vbic “ '^T ^ “ep^i ^ 

where the last term in this expression is Introduced to 
aooount foy majority carrier effect In the base side The 



relationships (II 8) to (XI 11) helps us determining 


and as a functionsof voltage To solve this equation 

first assume « 0 and solve for x^ With the value of 
P n 

X so obtained the same set of equation is solved for x 
n hr 
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itAiiuciion or bP Parameters 


/' Read in data f oni' 
CV mcasurcmcDt oO 
s. EB junciion 


Read in data frotti 
CV measurement on 
. CB junction 


Initial se the Ma qua dt 
paromete s and model 
pa omele s atues 


E t act the to elements 
of Pe related to £B 
ju ct on 


Extract the fo r elements 
of elated to C6 
) net on 


Read In the set 
FE ^c)of data 
gt onstant VcE&at loiv 
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